General method for calculating the universal conductance of strongly correlated 

junctions of multiple quantum wires 
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We develop a method to extract the universal conductance of junctions of multiple quantum wires, 
a property of systems connected to reservoirs, from static ground-state computations in closed finite 
systems. The method is based on a key relationship, derived within the framework of bound- 
ary conformal field theory, between the conductance tensor and certain ground state correlation 
functions. Our results provide a systematic way of studying quantum transport in the presence 
of strong electron-electron interactions using efficient numerical techniques such as the standard 
time-independent density-matrix renormalization-group method. We give a step-by-step recipe for 
applying the method and present several tests and benchmarks. As an application of the method, 
we calculate the conductance of the M fixed point of a Y junction of Luttinger liquids for several 
values of the Luttinger parameter g and conjecture its functional dependence on g. 

PACS numbers: 



I. INTRODUCTION 

Advances in molecular electronics can extend the lim- 
its of device miniaturization to the atomic scales where 
entire electronic circuits are made with molecular build- 
ing blocks*^ Single molecule junctions connected to two 
macroscopic metallic leads have already been successfully 
fabricated,— and there are several proposals such as laying 
quantum wires on top of each other for making junctions 
of multiple quantum wires^ 

If we eventually manage to build entire electronic cir- 
cuits with molecular building blocks, a paramount goal 
in the field of molecular electronics, junctions of three or 
more quantum wires will inevitably be a key ingredient. 
These junctions are comprised of several quantum wires, 
i.e., quasi-one-dimcnsional (ID) metallic structures with 
atomic scale sizes, that are connected to one another by 
a given molecular structure as shown schematically in 
Fig. [TJ The structure and interactions at the junction 
depend on the particular system under study. What we 
mean by metallic in the above description of quantum 
wires is that they are capable of conducting electricity 
due to the presence of gapless excitations. A generic 
description for these ID quantum wires is based on the 
Tomonaga-Luttinger-liquid theory^— Structures involv- 
ing Luttinger-liquid quantum wires have been the sub- 
ject of numerous recent studies.— ~— Experimentally, such 
quantum wires are realized with carbon nanotubes or 
through the cleaved edge overgrowth technique in GaAs 
heterostructures.— ~— Electrical current running in the 
wires can pass through this molecular structure at the 
junction. 

At molecular length scales, quantum mechanics is im- 
portant and the system represented above must be mod- 
eled accordingly. The simplest theoretical description of 
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FIG. 1: (Color online) A schematic illustration of the generic 
system that is the subject of this paper. We have M quantum 
wires connected to a molecular system of certain structure 
and interactions. We have currents I a running in wires a = 
1 . . . M and voltages V a applied to the endpoints of the wires. 



systems such as this is based on the tight-binding model, 
with anti-commuting creation and annihilation operators 
introduced for different atomic sites. An effective Hamil- 
tonian can then be written in terms of these creation 
and annihilation operators and generically involves hop- 
ping terms c^c and density-density interaction terms nn 
with n = c^c. 

Suppose we wish to study a rather arbitrary junction, 
the structure and interactions of which are represented 
by a tight-binding Hamiltonian. A very basic question 
regarding this system is how it conducts electricity. Con- 
sider a system of M wires. In the presence of voltage bi- 
ases V a applied to the endpoint of wire a for a = 1 . . . M, 
current I a will flow along each wire a. By convention, a 
current I a is positive if it flows toward the junction and 
negative if it flows away from it as seen in Fig. [T] In gen- 
eral, i.e., at arbitrary biases and temperatures, this prob- 
lem is very complicated and the currents flowing in the 
quantum wires are nonuniversal functions of the temper- 
ature, the voltages V a , and the microscopic details of the 
system. In this paper, we are concerned with the linear- 
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response regime where universal behavior can emerge. 
We work at zero temperature and consider only the limit 
of infinitesimal biases. The currents in this regime will 
be a linear combination of the applied biases as seen in 
the following: 

I a = ^2G a pVp. (1) 

P 

This linear relationship defines the linear conductance 
tensor G a p which is the quantity of interest in this paper. 

One of the most important ideas of modern physics 
is the remarkable universality that emerges near critical 
points. It turns out that due to the criticality of the 
bulk of quantum wires, i.e., having divergent correlation 
lengths, a large degree of universality also emerges in 
the behavior of quantum junctions. The universality can 
be understood in the framework of the renormalization 
group (RG).— One can argue that in the limit of small 
biases and low temperatures, many of the microscopic de- 
tails of the junction are irrelevant in the RG sense, which 
means their contribution to conductance, and other phys- 
ical observables, decays to zero at large distances and low 
energies. 

The junctions we are concerned with in this paper fall 
into the category of quantum impurity problems. The 
junction, with all the complex structure and interactions 
it contains, is localized at the endpoints of the wires. It 
can therefore be thought of as one (rather arbitrary) im- 
purity inserted into a system, the bulk behavior of which 
is given by that of M independent quantum wires. A clas- 
sic example of quantum impurity problems is the Kondo 
model describing the behavior of conduction electrons 
interacting with a local magnetic moment^ The power- 
ful methods of boundary conformal field theory (BCFT) 
have proven useful in a multitude of quantum impurity 
problems . 35 ! 36 Thus, BCFT is the main analytical tech- 
nique used in this paper. 

Determining the conductance of quantum junctions 
in the presence of strong electron-electron interactions 
is a long-sought and challenging goal. The Landaucr- 
Biittiker's formalism, which is the method of choice in the 
calculation of quantum conductance, does not account 
for these interactions, which indeed play a key role in 
low dimensions. Functional renormalization-group meth- 
ods have been helpful in studying the interaction effects 
in the vicinity of the junction, but their applicability is 
also dependent upon the presence of large noninteracting 
leads' 

In recent years, efficient numerical methods, such as 
the density-matrix renormalization group (DMRG), 39 
have been developed for studying strongly correlated 
quasi-lD quantum problems. Since the quantum junc- 
tions described above can be thought of as quasi-lD (by 
folding all the wires to one side so they run parallel to one 
another), these numerical methods could potentially be 
efficient tools for computing the conductance of junctions 
with an arbitrary number of wires and in the presence of 
strong interactions. 



In fact, DMRG has already been applied to the study 
of quantum junctions^ However, when it comes to cal- 
culating the conductance of strongly correlated junc- 
tions, there are fundamental difficulties even when we are 
armed with powerful tools such as DMRG. One such dif- 
ficulty arises from the fact that conductance is a property 
of an open quantum system. We define the conductance 
in terms of the current passing through the system and 
the underlying assumption is that we have reservoirs that 
can act as sources and drains for electrons. To study con- 
ductance, we either need to model the reservoirs carefully 
or send them to infinity. The latter is a simpler and more 
elegant way of formally dealing with quantum transport, 
but has the downside that for a numerical calculation of 
conductance, we would need to model large enough sys- 
tems that faithfully approximate the semi-infinite ones. 

Another difficulty with calculating the linear conduc- 
tance is that, within the linear- response framework, con- 
ductance is formally related to dynamical correlation 
functions. It may then appear that one needs to use the 
much more computationally demanding time-dependent 
numerical methods such as time-dependent DMRG to 
calculate the conductance. 

For junctions of two quantum wires, time-dependent 
DMRG has already been used for conductance calcula- 
tions4i~— A brute force calculation with time-dependent 
methods in large systems is not, however, currently fea- 
sible for strongly correlated junctions of more than two 
quantum wires. 

It is the objective of this paper to make such calcula- 
tions possible with a combination of analytical and nu- 
merical techniques. More specifically, the main objective 
of this paper is to develop a formalism that would allow 
us to apply numerical methods such as time- independent 
DMRG and the related matrix product states to calculate 
the linear-response conductance of strongly-correlated 
junctions of an arbitrary number of quantum wires with 
rather generic structures and interactions in the junc- 
tion^ In this paper, we focus on the systems with spin- 
less electrons, but our method can also be extended to 
systems with spin- 1/2 electrons. 

One particular application for the formalism we seek 
to develop in this paper is the problem of the M fixed 
point in a Y junction of spinlcss Luttinger liquids. The 
existence of this nontrivial fixed point was conjectured 
many years ago, but its nature, and more specifically its 
conductance, had remained an open question in quantum 
impurity problemSi 45 ' 46 

The outline of this paper is as follows. In Sec. HH we 
summarize the main results of this paper. We present 
a key relationship [Eq. ([2])] between the junction con- 
ductance and certain static correlation functions in a fi- 
nite system. This relationship serves as the basis of the 
method developed here. We also present in Sec. [TT] a 
step-by-step recipe for applying the method in practice 
as well as a summary of the new results on the Y junc- 
tion, which were obtained with this method. In Sec. IIII1 
we present some explicit calculations in a noninteract- 
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ing lattice model that help motivate the derivation of 
Eq. ([2} and clarify the connection of the continuum re- 
sults to lattice calculations. The results derived are in 
fact some special cases, which can be obtained with el- 
ementary methods, of the general relation Eq. ([2]), the 
derivation of which requires the machinery of BCFT. In 
Sec. lIVl we briefly review the main analytical techniques, 
namely, bosonization and boundary conformal field the- 
ory, used in this paper and set up the notation. Sec- 
tion [V] is devoted to deriving Eq. @ in the BCFT frame- 
work. In Sec. I VII we discuss in detail the method pro- 
posed in this work for conductance calculations and clar- 
ify practical issues regarding a lattice-model implementa- 
tion. In Sec. IVII1 we present numerical benchmarks with 
DMRG for interacting systems and exact diagonalization 
for noninteracting systems to verify the correctness of the 
method. In Sec. IVIII1 we study a Y junction of quantum 
wires and obtain the previously unknown conductance of 
its M fixed point as a function of the Luttinger parameter 
g. Finally, we conclude in Sec. IIXI by outlining the out- 
look for future applications and the impact of the results 
obtained in this paper. Some of the results of this paper 
have been briefly reported in Ref. 42. 



II. MAIN RESULTS 

Wc developed a method to extract the universal lin- 
ear conductance G a p of quantum multi-wire junctions 
defined in Eq. ([T]) from a calculation of the ground state 
expectation values of operators involving currents and 
densities in an appropriately constructed finite system. 

At the core of the method lies an important general 
relationship, which we recently derived in Ref. |47| using 
the machinery of BCFT. The relationship is derived in 
Sec. [Viand simply states that 
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where x is the distance from the boundary on the left in 
a system of length £, with I — > oo and finite x/£, con- 
structed from the junction of interest and an appropriate 
mirror image placed on the right endpoints of the wires 
as seen in Fig.[5J Here Jr(x) (jf (x)) is the right-moving 
(left-moving) current on wire a ((3). Note that although 
Eq. ^ holds asymptotically (x — > oo), it can be used to 
extract the conductance G a p even with finite but large 
enough x. 

The relationship above is the key ingredient of our 
method for calculating the conductance. If we can com- 
pute the quantity (JJj(x) J^(x))qs, which is a ground- 
state expectation value in a finite system of length £, as 
a function of x, then we can multiply it by a universal 
function to get the left-hand side of Eq. © above. This 
quantity will then saturate to the universal conductance 
of the junction for large x. 

Apart from the derivation of the key Eq. ([5]), we pro- 
vide a recipe for applying this continuum result to a lat- 
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FIG. 2: (Color online) A schematic illustration of the fi- 
nite system, constructed from the junction of interest and an 
appropriate mirror image, which we use in our method to ex- 
tract the conductance G a p. The conductance is related to 
the ground state ( Jr(x) jf (x))gs correlation function of chi- 
ral currents in this finite system through Eq. ([2]). 



tice calculation. This requires specifying the procedure 
for constructing the lattice (tight-binding) Hamiltonian 
of the aforementioned finite system from the Hamilto- 
nian of the junction of interest (that couples infinitely 
long wires). It also requires specifying lattice operators, 
i.e., in terms of the tight-binding fermionic creation and 
annihilation operators, the correlation functions of which 
are a good approximation to ( Jr(x) J^(x)}qs- This is im- 
portant because the chiral current operators are defined 
for the continuum theory, and chiral creation and anni- 
hilation operators can not be directly modeled on the 
lattice. The recipe for applying the key relation Eq. J2]) 
to a lattice computation is given in Sec. IVII 

For quick reference and an illustration of the method, 
here we give a simple example and explain a step-by-step 
application of the method to the well-known problem of 
a weak link in the Luttinger liquid. The starting point 
for applying our method is always a tight-binding lattice 
Hamiltonian of spinless electrons for bulk wires and their 
connection at the junction, namely, H = -Hboundary + 
-Hbuik- For a weak link in a Luttinger liquid, we can 
write 

-^boundary = —tc\ Q C2fi — tc 2 fiClfi 

and the following Hamiltonian for the bulk of the wires: 

2 oo 

ffbuik = [ ~ c «,j cq j+ i ~ c 1j+i c «j ( 3 ) 

a =l j=Q 

+V(n atj - -){n a , j+ i - i)]. 

Note that there is some arbitrariness in dividing the sys- 
tem into the junction and wires. The boundary Hamilto- 
nian above is a minimal choice, but including more sites 
in ^boundary does not affect the results as long as the 
system is large enough and the correlation functions dis- 
cussed below are computed far away from the boundary. 
Given the system Hamiltonian, our method consists of 
the following steps. 

1 . Construct the finite system as in Fig. [3J For this 
we need to construct a Hamiltonian H' = Hq + 
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H^ ulk + Hg where H and Hi respectively describe 
the junction on the left side of the system {x = 0) 
in Fig. [2] and the mirror image at x = £ . The recipe 
for constructing these Hamiltonians is simple. The 
left boundary Hamiltonian Hq is simply equal to 
-ffboundary and the bulk Hamiltonian -ff b uik nas ex ~ 
actly the same form as the bulk Hamiltonian of the 
semi-infinite system but a finite number of terms, 
i.e., Y^jLo ~^ Ylj=o ■ The construction of the right 
Hamiltonian goes as follows. First we consider the 
same Hamiltonian as ^boundary but acting on the 



other endpoint, namely — tc[ n C2,n — tct jy-c^jv and 
then we apply two transformations K and C on this 
Hamiltonian. K simply takes the complex conju- 
gate and C changes c — s- c) . In this case, assuming 
a real hopping amplitude t, we have 



Hi = C{-tc\ N c 2 ,N - tc\ N ci^ N ) 
= tc\ N c 2 ,N+tc\ N ci^ N . 

2. Having constructed the Hamiltonian H = Hq 



Hi 



bulk 



Hi of the finite system, measure, by a nu- 
merical DMRG calculation, the following ground 
state expectation value: 

(4(m)Jl(m)> = - 1 <ji( m )j2( m))) 



L a,m+1 



IS 



where J Q (m) = i ( c' a>m c a>m+1 - c 
simply the current operator and v is the charge car- 
rier velocity for the Luttinger liquids described by 
-ffbuik- The above equation is valid for time-reversal 
symmetric systems like the example at hand. The 
general construction of the operator {J^m) J^m)) 
in terms of the lattice creation and annihilation op- 
erators is given in section IVII 

3. Fit the data for ( Jp(m) J|(m)) to the asymptotic 
functional form from Eq. ©, i.e., ( J^.(?n) j£(m)) oc 

[4 N sin (t^-ti)] 2 , and obtain G\ 2 from the overall 
coefficient. 

Note that if Gyi = 0, the fitting is tricky. These details 
will be discussed later. 

The other important results of this paper concern a 
concrete application of the method to a previously un- 
solved quantum impurity problem. These results are 
presented in Sec. IVIIII The main problem solved in that 
section by an explicit numerical DMRG calculation is 
determining the conductance of a nontrivial fixed point 
in an interacting Y junction of three quantum wires de- 
picted in Fig. [3J The existence of this fixed point known 
as the M fixed point was conjectured in Ref. |45||46| . but 
the properties of the fixed point including its conduc- 
tance remained unknown. In this paper, by combining 
the method developed in Sees. IVl and fVTl and numerical 
computations with time-independent DMRG, we calcu- 
late the universal conductance of this fixed point as a 



function of the Luttinger parameter g (a parameter that 
quantifies the strength of electron-electron interactions 
in the wires). Based on the numerical results, we claim 
that the conductance of the M fixed point depends on 
the Luttinger parameter as 



G 12 (g) =G 21 (g) 



2ff7 



2.9 + 37 — 3<77 h 



(4) 



where 7 w 0.42. The above expression is a universal re- 
sult which, for any nonvanishing hopping amplitude t, 
holds independently of the value of t. We also explicitly 
verify that the conductance exhibits universal behavior. 
Namely we find that at large length scales, the conduc- 
tance is independent of the hopping amplitude t at the 
junction. 



Hopping amplitude t 



Magnetic flux 




Luttinger parameter g 



FIG. 3: A schematic illustration of the Y junction numeri- 
cally studied in section IVIIII The M fixed point corresponds 
to the time reversal symmetric case <j> = 0, ir. 



We also present in Sec. I Villi a numerical verification of 
a theoretical prediction for the conductance of the chi- 
ral fixed point. The chiral fixed point is stable in the 
range 1 < g < 3 and in the presence of a time-reversal 
symmetry breaking flux <fi ^= 0, 7r. 



III. CONDUCTANCE AND CORRELATIONS IN 
A NONINTERACTING LATTICE MODEL 

There is a well-established framework, namely the 
Landauer-Biittiker formalism (see for example Ref. HI) , 
for calculating conductances of multi-wire quantum junc- 
tions in the absence of electron-electron interactions. 
This framework is applicable both in the continuum and 
to lattice models. The key quantities in this framework 
are the incoming and outgoing scattering states "Fin and 
^out, which for a junction of M wires can be represented 
by an M x I (2M x I) column vector for spinless (spinful) 
electrons. The scattering states are labeled by momen- 
tum quantum number k, and the effect of the junction is 
encoded in an M x M (2M x 2M) unitary scattering ma- 
trix for spinless (spinful) electrons. The scattering matrix 
relates the incoming and outgoing scattering states as 



* out(fc) = S(k)V in (fe). 



(5) 
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In the Landauer-Biittiker formalism, conductance is 
simply related to the elements of the scattering matrix 
S(kp) at the Fermi level. Because including spin is a 
simple extension of the method developed in this paper, 
we work with spinless electrons here. The conductance 
between two different wires G a p is then given by 



G a /3 = \S a p{kp)Y 



(6) 



Let us now consider the simplest junction of two lattice 
wires. Each wire has hopping amplitude set to unity in 
the bulk and the two wires are connected by a hopping 
amplitude t. More generically, one can consider a nonin- 
teracting junction of M wires that are coupled quadrat- 
ically as described below. Consider a lattice system of 
M wires with electron annihilation operators c a j, where 
a = 1, M is the wire index and j = — oo, is the 
site index on each wire. We can write the Hamiltonian 



in a compact form by defining ^ 



(c! 



3 C 2,j 



C 



t 



u 



E 

i=—oo 



+ *Jr*o, (7) 



where r is a Hermitian M x M matrix with diagonal ele- 
ments r QQ = fx a (endpoint chemical potentials) and off- 
diagonal elements T a p = T*^ a = t a p (hopping between 
endpoints). In the bulk of each wire, the nearest- neighbor 
hopping amplitude is set to unity. The scattering eigen- 
states are 



U')k = *5 ( A+elk3 + A ~ e ~ ikj ) i°>> 



(8) 



and 
(+(-) 



where \&t is a row vector of operators c a 

is a column vector of scattering amplitudes AZ { 1 for 
a = 1 . . . M, and the standard matrix multiplication 
convention applies. By plugging the above states in 
the Schrodinger equation H\ip)k = ^k\ip)k, we obtain 
€k = — 2cosfc and^ 



s(k) = -(r + e - ifc )- 1 (r + e ifc ) 



(9) 



and certain static correlation functions in a finite sys- 
tem. The purpose of this section is to illustrate this 
relationship in the very special case of a simple nonin- 
teracting system where exact calculations can be done 
with elementary methods. Consider the following finite 
system shown in Fig. 01 which consists of our junction 
with hopping t on the left. At a finite length away from 
the junction, the other endpoints of the wires are coupled 
with a hopping amplitude — t. Let us now calculate the 
current-current correlation function (J 1 (a;) J 2 (x)). 




FIG. 4: A simple noninteracting system of two wires con- 
structed with a junction with hopping t and mirror image 
with hopping — t. The quantity of interest is the expectation 
value (J 1 (x)J 2 (x)) with J a the current operator on wire a. 

To begin with, let us assume a simple special case 
where t = 1. In this case, our system is a loop, with 
no impurity, that is threaded with a flux n. This gives 
rise to anti-periodic BC. In this special case, the Hamil- 
tonian can be simply diagonalizcd and we get 



H = — 2cos/c c^Cfc, k 



2n - 1 
2N ' 



where is the Fourier transform of Cj with j = 1, 



(12) 



.2N 



and n = -N + 1,...N. 
and assume N is even. 



We focus on the half-filled case 
The ground state is then given 



by a Slater determinant |GS) = c kl . ..c' k |0) and the 
fermionic correlation function 

C(a,b) = (GS|4 Cf) |GS), 

where a and b are the site index as in Fig. [5l can be 
written as 



The scattering matrix relates the incoming and outgoing 
states * in = A + and * out = A~ as * out = S^ hl . For 
the simple case of just two wires connected with hopping 
amplitude t, we have 



r = 



-t 

-t 



S 



t(e 2ik -l) 

e zk f2 _ £ — zk 

t 2 -i 

t 2_ e -2ik 



(10) 

The above scattering matrix at half-filling then yields the 
following conductance : 



-it 2 



(ii) 



As discussed in the Introduction, the key result of this 
paper is a general relationship between the conductance 



n=l 



c.c. 



(13) 



Each term in the above sum corresponds to one filled 
momentum level where the contributions of k and — k 
levels are complex conjugates. The sum is just a geomet- 
ric series, which can be calculated exactly and the result 
is 



C(a, b) 



1 sin 



7r(a— b) 



(14) 



sm ■ 



2N 



If J 1 (m) (on wire 1) is the current operator between 
sites m and m + 1, J 2 (m) will be between sites 2A^ — m 
and 2N — m+1 in the chain. Let us write these operators 
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explicitly as 

J 1 (m) = i (c^ C m+ i - 4„ +1 Cm) , 

</ 2 (to) = * (4w+l-m C2JV-m - 4jV-m c 2JV-m+l) • 

Using the fermionic Green's function Eq. ([Ft]) and Wick's 
theorem, we can explicitly calculate the current-current 
correlation function (J 1 (m) J 2 (m)). The real-space form 
of Wick's theorem can be generically written in the form 

( C ro C li °i C j) = ( C rn c j) ( c n c i) ~ ( C m c i)( c ri c i)- (15) 

To calculate (J^to)./ 2 ^)), we write it as a sum of 
I 



four quartic terms and reduce each term to a sum of 
products of single-electron Green's functions. Notice that 
C(a,b) = C(a — b), a function of the distance a — b alone. 
After some algebra, we can then write the current-current 
correlation function as 

{^{mj^im)) = 2C 2 (2N-2m) (16) 
- 2C(2TV - 2to + 1)C(2TV - 2?n - 1). 

By plugging the explicit form of C(a — b) from Eq. (|T4|) 
into Eq. (fT6"| , we then obtain the following exact expres- 
sion for the current-current correlation function: 



\j\m)J 2 {m)) 



snr 



[ir{N - to)] 



[ir(N ■ 



2N 2 



sin 2 [f (TV -to)] sin [f (N ■ 



2N J 



[tt(TV- 



2N J 



Consider the above correlation function away from the 
two endpoints TV — to 3> 1 and to 3> 1. We can 
then approximate the denominator of the second term 



as sin 2 [jf(N — to)] and write 



(J^to^J^to)) 



1 r 



TV sin 



/ to \ 

("tv) 



(17) 



As we will show later, the expression derived above for 
a simple nonintcracting model has a universal form that 
survives electron-electron interactions. 

Before proceeding, let us consider a limit of the above 
expression. If we take the limit of TV — > oo in Eq. (fT7|) 
above, we are sending off the right junction to infinity 
and effectively describing a semi-infinite system. In this 
case, for a distance to (lattice spacing set to unity), we 
obtain the following correlation function: 



(J 1 (m)J 2 (m)) = -LJ^. 

Z7T Z TO 



(18) 



In the next step, we derive a similar expression for a 
semi-infinite system but with an impurity consisting of a 
site with hopping amplitude t. It is illuminating to first 
give a short derivation of Eq. (|18l) formulated directly in 
the limit of the semi-infinite system. In this limit, we can 
treat the momenta in the continuum and work directly 
with the scattering wave functions, which in the special 
case of t — 1 are just plane waves e lkx due to the absence 
of back-scattering. We have right-moving (left-moving) 
plane waves for < k (0 > k) and at half-filling, the filled 
momentum states have — ? < k < This leads to the 
following fcrmion Green's function: 



ika ^—ikb 



7r(a— b) 



ir(a — b) 



(19) 



which is in fact the limit of TV — > oo of Eq. (|14p . Inserting 
the above expression into 

(J^to)^^)) = 2C 2 {a-b)-2C(a-b+l)C{a-b-l), 

where m — (b — a)/2 straightforwardly leads to Eq. (|18p . 

So far, we have only considered t = 1. Let us now 
consider a semi-infinite system shown in Fig. [5] with an 
impurity of arbitrary hopping amplitude t. Now, instead 
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FIG. 5: (Color online) The single-particle scattering states 
for a single-impurity infinitely long system. 

of the simple plane waves e lkx , we can use the following 
single-particle scattering states. Let us write the right- 
moving and left-moving scattering states separately for 
clarity. We have 



4$(x) = 



e ikx _|_ rfe g 

t k e ikx , 



- lkx , x sC 
x > 1 



e lkx + f k e lkx , 



the 



x > 1. 



In terms of the above scattering states, we can write 
the fermionic Green's function for a < and b > 1 as 



C{a,b) 



+tle 



* dk 
2tt 

ika 



[{e~ lka + r* k e tka ) t k e 



ikb 



-ikb 



f k e ikb 



)]■ 



(20) 
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Note that in the simple case t = 1, where we have no 
back-scattering and r k = r k = 0, Eq. (f2TJ)) simply re- 
duces to Eq. (|T9l) above. By using the scattering matrix 
Eq. (|10[) or more directly from plugging in the scattering 
states into the equations 

e*^fc(l) = -*^fc(0) - 0*(2), (21) 
e fc <M0) = -t&(l) - &(-!), (22) 



we obtain = 
< fe = t k and f fe 



-2ifc_ t 2 

-2ik 



and t k =t- 



We also have 



-2ik_ t 2 ■ 

e ~""r k . By inserting the above trans- 
mission and reflection coefficients into the expression for 
the Green's function Eq. (f2U|) and after some algebra, we 
can write 



C(a, b) = At 



2 dk , sin \k(a — b 

— smfc 

2tt 1 



l)]t 2 -sm[k{a-b-l)} _ 



2t 2 cos 2k + t 4 



= C(a-b). 



r 



(23) 



This leads to the following expression for the current- 
current correlation function: 

(^(mj^im)) = 2C 2 (-2m) - 2C(-2m - l)C(-2m + 1) 

(24) 

with C(x) defined in Eq. ([23]) . We claim that the asymp- 
totic behavior of the quantity {J 1 {m) J 2 (m)) is generi- 
cally given by ~ -K as in the special case Eq. (fl"8|) . To 
check this claim we plot m 2 (J 1 (m) J 2 (m)) as a function 
of m for several values of the hopping amplitude t by 
straightforward numerical evaluation of the integral in 
Eq. (|23|) . The results are shown in Fig. HI 



behavior: 

(J l {m)J 2 (m)) 



At 2 



2tt 2 (l + t 2 Y 



m > 1. (25) 



At this level, this observation appears rather myste- 
rious. It is not very transparent from the expression 
for (J 1 (m) J 2 (m)), which is a rather complicated dou- 
ble integral, why this ground state correlation function 
has such a simple scaling form. By comparing Eq. (|25[) 
above with Eq. (|17[) . we claim that for a finite system 
shown in Fig. 21 the asymptotic behavior of the correla- 
tion function is given by 




FIG. 6: (Color online) The quantity m 2 (J 1 (m) J 2 (m)) calcu- 
lated from Eq. (|24[) by a numerical integration of the integral 
in Eq. I|23[). For large m, this quantity saturates to a constant 
value that only depends on the hopping amplitude t. This in- 
dicates that the asymptotic form of (^(m) J 2 (m)) is given by 



2V 



(i + t 2 y 



iVsin 



m 7T — 

V NJ 



m > 1. 



(26) 

We do not attempt to prove the above expression here. 
The equation is indeed a very special case of the generic 
relationship Eq. ([2} that we prove in section |V] in full 
generality, i.e., in the presence of electron-electron inter- 
actions for a rather arbitrary junction with an arbitrary 
number of wires and without assuming symmetries such 
as time-reversal. 

The expression above, however, motivates the main 
idea of this paper, i.e., the fact that conductance can 
be extracted from ground state expectation values in a 
closed system, in a very elementary example. It is also 
an example where the relationship can be derived explic- 
itly on the lattice instead of resorting to the continuum 
formalism, which helps clarify the application of the con- 
tinuum CFT results to a lattice calculation. 



Since for each t, m 2 (J 1 (m) J 2 (m)) saturates to a con- 
stant, a natural question is how this constant depends on 
the value of t. A quite remarkable fact is that the sat- 
uration value is exactly proportional to the conductance 
Eq. (jlip of this simple junction, which we calculated in 
the beginning of this section from the Landauer's formal- 
ism. In other words, we have the following asymptotic 



IV. REVIEW OF BOSONIZATION AND 
BOUNDARY CONFORMAL FIELD THEORY 

The main results of this paper are derived within the 
framework of boundary conformal field theory. In this 
section, we review the required steps for formulating the 
generic system we would like to study, namely a junction 
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of M quantum wires modeled as a tight-binding Hamil- 
tonian with spinless electrons, in the language of the con- 
formal field theory. The material in this section can be 
skipped by readers familiar with bosonization and CFT. 



A. Bosonization of quantum wires 

The first step in this formulation is the bosonization 
procedure. Bosonization is a powerful nonlocal transfor- 
mation that, in one space dimension, allows us to describe 
the low-energy limit of strongly-correlated fcrmionic sys- 
tems as a noninteracting theory of bosons (see, for exam- 
ple, Ref. HI for a review) . Let us start by considering one 
infinitely long quantum wire with the following Hamilto- 



-c\ci + i - c} 



+ V(th- - i) 



(27) 

where n, = c|cj. The first two terms describe electron 
hopping and the last term is the density-density inter- 
action. At half-filling, this model exhibits a charge den- 
sity wave phase transition for large repulsive interactions 
V > 2. We then see that the ground state spontaneously 
breaks lattice translation symmetry and we get the two 
degenerate ground states shown in Fig. [7] for V — > oo. 



FIG. 7: The degenerate ground states of the Hamiltonian 
Eq. (|27[) for V — > oo. Filled circles represent occupied sites. 

Also at very large attractive interactions, the electrons 
will clump together and form clusters of neighboring oc- 
cupied sites. This results in phase separation which, at 
half-filling, happens at V < —2. 

For a certain range of interactions — 2 < V < 2, which 
includes the noninteracting fermionic system (V = 0), 
the system is in a gapless critical phase known as the 
Luttinger liquid. This is the regime where bosonization 
applies. Let us first bosonizc the noninteracting system. 

The noninteracting fermions have an tk = — 2cos(/c) 
dispersion and a ground state consisting of a filled Fermi 
sea up to the Fermi level hp. We can linearize the disper- 
sion around the Fermi level and only consider the right- 
moving (left-moving) excitations close to k F (—k F ). 

Let us define the following right-moving and left- 
moving fields: 

ip R (p) =c(k F +p), ip L (p)=c(-k F +p). (28) 

We can then write the hopping part of the Hamiltonian 

as 



Ho 



dp 

VF 1 Y 



P>Pr(p)^r(p) - PipLipjipLip) • (29) 



Upon Fourier transforming, the above equation in real 
space reads as 



Ha = vp dx 



where 



Id Id 



ik F x 



(30) 



(31) 



The key step to bosonization is writing the Hamilto- 
nian (both hopping and interaction part) in terms of the 
following chiral current operators: 



Jr,l(x) = ij) R L {x) ^r,l(x) 



(32) 



instead of the fermionic creation and annihilation opera- 
tors. One can easily verify using the Hamiltonian Hq in 
Eq. ([3H)) that these currents satisfy the following commu- 
tation relations: 



[H ,j Rt . 



+ • 9 ■ 
±iv F —j RtL . 

ox 



(33) 



The interaction term, which is quartic in creation and 
annihilation operators, will trivially become quadratic 
in terms of the currents. The nontrivial part of the 
bosonization procedure is to show that the hopping part 
of the Hamiltonian, which is quadratic in -0, will remain 
quadratic in terms of j. 

This important result can be shown using the com- 
mutation relation of the current operators. By Fourier 
transforming tpR as ^Pr(x) = e V'B(fc), we can 

write the Fourier transform of the chiral current jr in 
Eq. fl32J) as 



3r{i) 



q)ipR{k), 



(34) 



which leads to the following commutation relation^ 



50. 



\jn(q),jR(<l')] 



q A 



(35) 



A similar equation with j- — > — can be derived for . 
Using the commutation relation Eq. Q35p. one can show 
that the bosonized Hamiltonian 

H = nvp \jR(q)jR(-q) + 3l{i)3l{-i)\ (36) 



obeys the commutation relations in Eq. (f3"3"|). We then 
argue that Hamiltonian (|36p is the bosonized form of the 
noninteracting fermionic Hamiltonian (|30[) . 

By using the commutation relations between j_r.l and 
i>R,L , one can write the fcrmionic operators in terms 
of the currents by introducing bosonic fields §r,l- The 
result is 



1 



2tt 



3r,l = ^-S-(pR,L- (37) 

27T OX 
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The interacting Hamiltonian density will now have di- 
agonal terms of the form jf^(x), j\(x) coming from the 
hopping and off-diagonal terms of the form 3L{x)jn{x) 
from the electron-electron interaction (back-scattering) 
and is quadratic in terms of the currents. By perfroming 
a linear transformation 



Jr 
Jl 



cosh /3 sinh (3 
sinh (3 cosh /3 



(38) 



that preserves the commutation relations, we obtain 

U{x)<x[J 2 R {x) + J 2 {x)]. (39) 

It is convenient to define new bosonic fields tp and 9 
that are linear combinations of 4>r,l in Eq. (|37[) such that 



iPr,l( x ) = 



1 



J( V ,(x)±S(x))/V2 



(40) 



In terms of these fields and by introducing the renormal- 
izcd charge carrier velocity v and the dimcnsionlcss Lut- 
tinger parameter g, the low energy effective Hamiltonian 
density for a wire can be written as 



U(x) 



v 

4tt 



f 



g(d x ipY + -(d x 6) 
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(41) 



The bosonized Hamiltonian Eq. (|4"1"T) is the generic de- 
scription of ID metallic quantum wires that we work with 
in this paper. Notice that ^-d x tp is the momentum con- 
jugate to 6 and we have the following commutation rela- 
tions: 



[<f{x), 0(x')] = iwsgn(x' — x). 



(42) 



The above review covers the main ingredients of the 
bosonization scheme that we need for setting up the prob- 
lem at hand in this paper. 



B. Boundary conformal field theory 

Here, we briefly review the basics of CFT and BCFT 
that we need in the remainder of this paper. This review 
is largely based on Ref. l5lT An important property of 
critical theories, i.e., theories where certain fields known 
as the scaling fields have critical correlation functions, is 
scale invariance. This simply means that the correlation 
functions of scaling fields 0\ satisfy 



(0 1 (br 1 )0 2 (br 2 )...0„(br n )) 

= b-^ x *(0 1 (r 1 )0 2 (r 2 )...O n (r n )), 



(43) 



where the exponent x\ is known as the scaling dimension 
of the operator Oi and b is an arbitrary scaling factor. 

A powerful leap from scaling symmetry to conformal 
symmetry is by allowing the scaling factor b to vary 
smoothly, i.e., considering more general transformation 
r — >• r' than r — > br with b(r) = |t£t|. If we have a local 



theory and the transformation r — > r' locally resembles 
a scaling transformation (modulo a local rotation) then 
we expect 

(O l (r[)0 2 (r' 2 )...O n (r' n )) 

= II b ( r ^~ Xi ( 0l (n)°2(r 2 ) . . . O n (r„)>. (44) 

i 

In three dimensions, the transformations that locally re- 
semble scaling are limited but in two dimensions (2D) 
(or 1+1D quantum systems) any analytic mapping [z — > 
w(z) on the complex plane z] is conformal (preserves an- 
gles) and requiring conformal symmetry leads to highly 
nontrivial results. An important quantity in a CFT is the 
stress-energy tensor. Consider the action S of a 1+1D 
quantum system that is defined on the complex plane z. 
The stress-energy tensor is defined in terms of the 
variations of this action through 



SS= -— I T^d v a»d z r, 



(45) 



where 8S is the variation of S due to the infinitesimal 
transformation r M — > r M + a^(r). Another important 
ingredient of CFT is the operator product expansion 
(OPE), which describes the nature of the singularities 
in {O l {r i )O j (r j ) . . .) as n ->• r 3 as 

k 

(46) 

In any CFT, the primary fields are fields for which the 
most singular term in the OPE of T(z)Oj(Zj, Zj) is order 
(z - Zj )- 2 where T(z) = T zz . 

Let us consider primary fields with the following scal- 
ing transformation on the complex plane z: 

Oj{\z,\z) = \- A t\- K iOj{z,z), 

where Aj and Aj are called the complex scaling dimen- 
sions. Notice that the above equation is only a shorthand 
to describe the scaling behavior of correlation functions 
involving O. An important result in CFT is that under 
a generic conformal transformation z — > w(z), the corre- 
lation functions of primary fields transform as^i 



(Oi{zx, z 1 )0 2 (z 2 , z 2 ) 

=n 



duii 


Ai 


dwi 


dzi 




d^ 



A, 



(Oi(wi,«;i)...)- (47) 



This is in fact the same as Eq. (|44[) which we intuitively 
wrote down in the beginning of this section. We shall 
emphasize that the transformation (|4"T|) only holds for 
the primary fields of the theory. 

So far our discussion has been limited to CFTs on 
the entire complex plane z. A CFT on a domain with 
boundaries is known as boundary conformal field the- 
ory (BCFT) where in addition to the bulk properties, we 
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need to specify appropriate boundary conditions (BCs). 
An important task in BCFT is classifying the conformally 
invariant boundary conditions for a given bulk CFT.— 

The BCFT techniques have proved powerful in study- 
ing various quantum impurity problems such as the 
Kondo model and Luttinger liquids with impurities ! 36 ! 53 
The effect of the impurity, in this approach to quantum 
impurity problems, is to select the appropriate confor- 
mally invariant BC corresponding to the low-energy RG 
fixed point. At far away from the boundaries, the bulk 
CFT is expected to adequately describe the system, i.e., 
the boundary conditions do not matter. Also very close 
to the boundary, there is nonunivcrsal short-distance 
physics from the microscopic degrees of freedom. On 
distances from the boundary that are much larger than 
the microscopic length scales, but still much smaller than 
the domain size, the correlation functions are affected by 
the conformally invariant boundary conditions^ 

As pointed out by Cardy, the BCs in a BCFT are en- 
coded in boundary states^ To understand the notion of 
a boundary state, it is helpful to consider the partition 
function of a CFT on a cylinder with boundary condi- 
tions A and B on two sides of the cylinder as in Fig. [5^ 
One can write down the partition function using transfer 
matrices running parallel or perpendicular to the bound- 
ary. The transfer matrix in the direction parallel to the 




(' 



FIG. 8: Boundary states in the partition function of a CFT 
on a cylinder. 

boundary depends on the boundary conditions A and B 
and can be written as exp (—Hab) (the lattice spacing 
is set to unity). The imaginary time runs parallel to the 
boundary and H-ab is the Hamiltonian with boundary 
conditions A and B. We then have the partition function 

Z = tr [e-? HAB ] . (48) 

Alternatively one can construct the transfer matrix in the 
direction perpendicular to the boundary. The Hamilto- 
nian H for this transfer matrix is determined by the bulk 
only and is independent of the boundary conditions. To 
write the partition function we need to construct states 
\A) and \B) in the space where the transfer matrix acts 
such that the partition function Eq. (|4"5|) is also equal to 

Z = (A\e- £H \B). (49) 

The two different representations of the partition func- 
tion above illustrate the notion of the boundary stated 



V. A GENERAL RELATIONSHIP BETWEEN 
THE CONDUCTANCE AND STATIC 
CORRELATION FUNCTIONS 

Let us recall the two difficulties with calculating the 
conductance of a junction of multiple interacting quan- 
tum wires, which we discussed in the Introduction. First, 
conductance is a property of an open quantum system. 
This means that we think of the wires emanating from the 
junction as being attached to reservoirs that can serve as 
sources and drains for the electrons passing through the 
molecular structure at the junction. To study the con- 
ductance, we either need to faithfully model the reser- 
voirs in our theoretical description or assume that we 
have infinitely long wires. In the latter case, we basically 
send the reservoirs to infinity. An appropriate boundary 
condition is assumed at infinity, which does not enter the 
theory of the junction explicitly. For a numerical calcula- 
tion of conductance, it then seems that we would need to 
model large enough systems that approximate the semi- 
infinite ones. 

The second difficulty is that, in the standard linear- 
response framework, the linear conductance is related to 
dynamical rather than static correlation functions. This 
can be seen explicitly from the Kubo formula for the 
conductance tensor— 

e 2 1 f°° 

G a0 = lim - — / dre lu)T (50) 

x f dx (TrJ a (y,T)J?{x,0)), 
Jo 

where T T indicates imaginary-time ordering and the 
quantity (T T J a (y,T)j"(x,0)) is a dynamical current- 
current correlation function for currents J a and J* 3 on 
wires a and /3, respectively. It seems that, gencrically, 
a numerical calculation of this quantity requires time- 
dependent methods such, as for example, time-dependent 
DMRG. 

With a brute force approach, the numerical calcula- 
tion of conductance is not feasible for junctions of three 
or more wires. In this section, we derive a key relation- 
ship between the universal conductance and certain static 
correlation functions in a finite system. This important 
relationship is the basis of a method for numerically cal- 
culating the conductance that we will discuss in Sec. I VII 



A. CFT for independent wires 

Let us begin by formulating the problem in the BCFT 
framework reviewed in Sec. [IVJ We have M quantum 
wires connected to a junction. Let us first consider the 
system of M independent quantum wires in the absence 
of the junction. As argued in Sec. IIV Al the effective 
Luttingcr-liquid Hamiltonian for one wire can be written 
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H 



-. — / dx 

4ir 



~(8 X 6) 2 
9 



(51) 



Using the fact that lie = j^d x tp is the conjugate mo- 
mentum of 9, the Lagrangian density can be written as 



C 



4irg 



-{Of -v(8 x 6) 2 



Let us set the charge carrier velocity v to unity for con- 
venience. We can then write the Euclidean action as 



S = I drdxdjd 1 
4ng 1 



(52) 



Alternatively the action can be written in terms of the 
dual field <p as S — j- J drdx d^ipd^ip. It is convenient 
to represent the points on the [x, r) plane with complex 
coordinate z = r + ix. The system described by the 
massless action Eq. (|52[) is a CFT on the complex plane 
z. Physically, covering the whole complex plane implies 
ID quantum wires extending from — oo to oo at zero tem- 
perature (/3 — > oo). 

Let us begin by calculating the correlation functions of 
the bosonic fields 9 from the action Eq. (f5"2"]). We have 



(%9(X U T 1 )9(X 2 ,T 2 )) 



J V[9] 9(x 1 ,T 1 )9(x 2 ,T 2 )e- s 

JV[9] e-s 
K(xi—X2,n—n), (53) 



where the propagator K(x\ — x 2 , T\ — r 2 ) is the inverse 
of the operator — ^ + i.e., 



2wg 

and is given by 



(flg + 8 2 ) K(x -x',t- t') = 5{x - x')S(r - r') 

(54) 



K(x-x',t-t') = - -In [(x- x') 2 + (t-t') 2 ] + const. 

(55) 

Up to an unimportant additive constant, we can then 
write the correlation function of the bosonic fields 9 as 

(%9(zi,zi)9(z 2 ,z 2 )) = -- [ln(zi - z 2 ) + ln(zi - z 2 ) ] . 

(56) 

By differentiating the above Eq. (|56p with respect to z 
and z, we obtain the following correlation functions: 



^ 89 , _ , 86 , _ s 
It -^—{zi,zx) -^—{z 2 ,z 2 ) 
8z\ az 2 



2 [zi - z 2 y 



^ 89 86 \ g 1 



8z\ 



8z 2 



Also we can similarly show that {7~ T -§§^(zi, zi)-^(z 2 , z 2 )) 
vanishes. 



2 (zi - z 2 ) 2 



.(58) 



The chiral currents below, where we have put in a wire 
index a, are then primary fields for this CFT: 

JZ(z) = 9 a (z, z), J«(f) = --L. 5 6 a (z, z), 

(59) 



with the notation 



8 = 8 Z = -(8 T 



id x ), 8 = 8 S 



id x )- 



In the absence of a boundary, the chiral currents in 
different wires arc uncorrelated and the only correlations 
are between chiral currents of the same chirality in the 
same wire. The only nonvanishing correlation functions 
of the chiral currents are then 



(% ji ( Zl ) 4 ( Z2 )) = * 7 _** 

47H [z\ - z 2 y 



Sat 



(60) 



4?T 2 (Z! - z 2 ) 2 ' 



Note that the fact that there is no correlation between 
the left-movers and the right-movers does not mean that 
there is no back-scattering in the system. Indeed, we 
actually have back-scattering in the presence of electron- 
electron interactions. The chiral currents Jl,r are the 
chiral eigenmodes of the interacting theory which, as ex- 
plained in Eq. (|3"5)) , arc in fact linear combinations of the 
bare (nonintcracting) chiral currents jh.R- Interactions 
induce back-scattering for these bare currents jl.r- 

The form of the correlation functions above is highly 
constrained by conformal symmetry and is basically de- 
termined by the scaling dimension of the primary opera- 
tors. A very important observation is that space and time 
are tied together in the above correlation functions and 
a measurement of static correlation functions uniquely 
determined the dynamical ones. Let us note that these 
chiral currents are related to the physical current J and 
the density fluctuation p through 



J° 



- T° 



J? 



The velocity of the charge carriers v coincides with the 
Fermi velocity only in the absence of electron-electron 
interactions. For simplicity, we shall set v to unity in the 
remainder of this section. 



B. BCFT for the semi-infinite system 

Let us now return to the problem of the quantum junc- 
tion. The system of M independent wires is described by 
a CFT of M bosonic fields living on the full complex z 
plane. The junction of the M quantum wires at the RG 
fixed point is described by BCFT on the upper-half com- 
plex plane. The real axis x = is the boundary of the 
domain where the CFT lives. We expect to have the 
same bulk CFT as the system of M independent wires if 
we are infinitely far away from this boundary. 
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This is an example of a quantum impurity problem. An 
important hypothesis, which has been repeatedly verified 
in a multitude of quantum impurity problems,— such as 
the single-channel and multichannel Kondo model^ is 
that at the RG fixed points, the conformal symmetry in 
the bulk terminates smoothly. More precisely, the hy- 
pothesis states that the BC on the boundary of the do- 
main in the complex plane is conformally invariant and 
therefore the quantum impurity problem at the RG fixed 
point is described by a BCFT. 

How can the presence of the boundary change the cor- 
relation functions? In other words, what is the difference 
between the correlation functions in the CFT on the full 
z plane and the BCFT on the upper half plane? A key 
observation is that the presence of the boundary does not 
change the correlation functions between chiral currents 
of the same chirality. On physical grounds, it is easy to 
see that the left movers coming from infinity toward the 
junction have not yet felt the presence of the junction 
and, just by causality, can not have different correlation 
functions than in the system of M independent wires. For 
the right movers, on the other hand, we need the confor- 
mal invariance of the boundary conditions to make such 
a statement . 

If we have a conformally invariant BC on the real axis, 
(Jr(z~i) Jr(z 2 )) will be the same as on the full complex 
plane. This is generically the case in the limit when the 
two points Z\ and Z2 are far away from the boundary and 
deep in the bulk of the system. With conformally invari- 
ant BCs, however, these points can be brought close to 
the boundary through conformal transformations, such 
as z — > Xz or z — > — l/z, from the upper half-plane onto 
itself, which leaves (Jr(zi)Jr(z 2 )) invariant. The above 
argument is not a rigorous proof, but may give some intu- 
ition. Notice that irrelevant boundary operators can give 
corrections to (Jr(zi)Jr(z 2 )) m the semi- infinite system, 
which decay as a power law with the distance from the 
boundary with exponents larger than 2. In a finite sys- 
tem, both (J£(zi)J£(z 2 )) and ( J%(zi) Jr{z 2 )) will have 
such corrections. 

The presence of a boundary, however, does introduce 
new correlations between the left and the right movers. 
From the scattering picture, we can understand this in 
the sense that left movers must go through the junction 
before turning into right movers. The form of the new 
correlation functions is determined by the scaling dimen- 
sion of the current operators, and as far as these corre- 
lation functions are concerned, all the information about 
the boundary condition at x = is encoded in an overall 
coefficient A^f . 

In summary, in the BCFT describing the junction at 
the RG fixed point, in addition to the correlation func- 
tions Eq. (|60[) . we have the following nonvanishing corre- 
lation functions^: 



and a similar equation for (J~ T (z\) Jr(z 2 )). 



(61) 



The coefficients A^ are determined by the confor- 
mally invariant boundary condition on the real axis. In 
fact we can write this coefficient in terms of the boundary 
state \B) as££££ 



4^ B 



(0\B) 



(62) 



where \0, 0) is the highest weight state corresponding to 
a generic operator O (here O = J^Jr) an d |0) is the 
ground state. The definition of the highest weight state 
and a derivation of Eq. (|62p is given in appendix [B] 

One important observation is in order. While the 
Hamiltonian is local in terms of the fermionic degrees 
of freedom, the mapping from fcrmion to bosons, i.e, 
the bosonization procedure, is highly nonlocal. One 
must notice, however, that regardless of the number of 
wires, there is a fundamental difference between a two- 
dimensional system and this quasi-lD system. The non- 
locality only shows up along the boundary itself. The 
statistics of the microscopic degrees of freedom can af- 
fect what the boundary condition B, and consequently 
the coefficient A°g is but it does not change the fact 
that the effect of the junction can be reduced to a BC at 
x = 0. The generic form Eq. (|5T|) holds even though the 
bosonic fields arc nonlocal objects in terms of the original 
fermions. 

So far, we have established that the current-current 
correlation functions appearing in Eq. (|50[) have a univer- 
sal form modulo coefficients Ag 13 , which depend on the 
boundary conditions. We should then be able to express 
the conductance G a p in terms of these coefficients A°g 
by performing the integrals in the Kubo formula. Con- 
sider the conductance G a p between two distinct wires 
a ^ ft. Using J = Jr — Jl (velocity set to unity), we can 
write 



(T T J a {z 1 ,z 1 )J fi {z 2 ,z 2 )) 



y 

4tt 2 



.4 



a/3 



{zx-z 2 f 



(zi - z 2 ) 



Note that the presence of the boundary does not 
change the correlation between currents of the same chi- 
rality, and the chiral correlation functions are obviously 
equal to zero in the half-plane for a ^ (3 (because they 
are proportional to 5 a p). Plugging the above equation 
into the Kubo formula Eq. (|50[) (with z\ = r + iy and 
z 2 = ix) gives two similar terms. Let us first perform the 
integral over the imaginary time r. We can write 



dr 



(t — iu) 2 



-2moH(u)e 



ii^O (63) 



which can be easily derived by a contour integration. 
Here H(u) is the Hcaviside step function H{u) = 1(0) 
for u > 0(it < 0). Using the above identity we can then 



13 



write 

G a p 



lis m 



dx 



Af H(x + y) + A^ H(-x-y) 



which simply yields 



G 



a/3 



A af) 



9 ~h' a ^@- 



(64) 



In the above expression, since both x and y are posi- 
tive (the wire extends from x = to +oo), the second 
term identically vanishes. Notice that determining the 
off-diagonal elements G a p with a ^ (3 is enough to also 
determine the diagonal elements of the conductance ten- 
sor G. The elements are not independent and satisfy the 
relations 



^ G a /3 — 0, ^ G a (j — 0. 



(65) 



and, as a function of the distance x from the boundary, 
fit it to the power-law form above, we can extract the 
coefficient Ag 3 and uniquely determine the conductance. 

Because the BCFT description applies only at low en- 
ergies and large distances, the behavior of the system at 
short distances (in units of the lattice spacing) is nonuni- 
versal and depends on the microscopic details rather than 
the continuum limit. Even at distances where the con- 
tinuum limit is valid, the expression above only gives 
the asymptotic behavior of the correlation function. The 
conformally invariant boundary condition B describes the 
RG fixed point, but, generically, there are corrections to 
the BCFT prediction Eq. (|66|) from irrelevant boundary 
operators at finite distances. These corrections of course 
die out faster than , and at large distances the asymp- 
totic behavior is given by the leading term Eq. (|66[) . 



Conformal mapping from the semi-infinite to a 
finite system 



The first relation follows because applying the same volt- 
age to all the wires leads to zero current, and the second 
relation follows from current conservation. 

Notice that in deriving the above expression (|6"4"|) . we 
assumed that the interacting Luttingcr-liquid wires ex- 
tend to infinity. In many experimental situations, the 
interacting wires are attached to noninteracting Fermi 
liquid leads. If the interacting region is long enough, 
the effect of the Fermi-liquid leads can be taken into ac- 
count by considering a contact resistance as discussed 
later in Sec. IVIIII Short interacting wires attached to 
Fermi-liquid leads can be considered as part of the junc- 
tion, while leads constitute quantum wires with g = 1, 
which extend to infinity. The possible crossover behavior 
as a function of the length of the interacting region in the 
presence of infinite noninteracting leads is an interesting 
open question. 

One of the two difficulties with a numerical calcu- 
lation of the universal conductance is now effectively 
solved. The conductance is uniquely determined by the 
coefficient in front of the dynamical correlation function 
Eq. (|61[) . But, conformal symmetry determines the uni- 
versal form of this correlation function. In other words, 
space and time are tied together by conformal symme- 
try, and the coefficient A^f can be extracted from a 
static correlation function alone. More explicitly, con- 
sider the static correlation function {J^{x)J^{x)) that is 
merely a special case of the generic correlation function 
(T T Jft(z~i) J^{z-i)) in Eq. (|6T]) for z\ = z 2 = ix and is 
given 



(Jr(x)4(x)) 



9 A a/3 



At: 2 



1 



(2s) 



(66) 



The above correlation function is just a ground-state ex- 
pectation value in a semi-infinite system. If we can nu- 
merically compute this ground-state expectation value 



Having solved one of the difficulties of calculating the 
conductance, we now turn to the other difficulty. Namely, 
the fact that modeling an open, i.e., semi- infinite, sys- 
tem requires finding correlation functions in very large 
systems. Let us go back to the BCFT description of 
the semi-infinite junction for which the conductance is 
well defined. We argued that this system is described by 
a BCFT on the upper half-plane. As explained below, 
we can use a conformal mapping to go from the upper 
half-plane to a strip of width £, which describes a finite 
system of length £ at zero temperature. We will show in 
what follows that the correlation function (J^(x)J^(x)} 
in this finite system also has a universal form and allows 
the extraction of the key coefficient A^f . 

Let us first explain the mapping, shown in Fig. [51 from 
the upper half-plane z = r + ix to a strip w = u + iv. 
The mapping is given by 



w 



In z, 



(67) 



In polar coordinates we can write z = re 1 ^ on the half- 
plane with < r < oo and < <p < ir which gives 
u = — lnr and v = ^(j>. The x > part of the real axis 
in the z plane is mapped to the boundary v = and the 
x < part to v = £. 

In the finite system of length £, which we obtain from 
the conformal mapping above, the correlation functions 
behave in a universal manner. Let us formally consider 
the static correlation function of Eq. (|66| , which is only 
a function of the distance from the boundary x, as a 
function of z = t + ix but with no dependence on r. If 
we have a primary operator O (here O = J%J^ L ) in the 
semi-infinite plane with 



(0{z))=A°{2x)- x °, 



(68) 
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FIG. 9: (Color online) The conformal transformation from 
the upper half-plane to a strip of width i. The physical sys- 
tem described by the half-plane is a semi-infinite junction at 
zero temperature with the imaginary time running from — oo 
to +00 and the position along the wires running from to 
+00. The physical system corresponding to the strip is a fi- 
nite system of length £ at zero temperature. The structure of 
the junctions at the two end-points of this finite system must 
give rise to the appropriate boundary conditions consistent 
with the conformal mapping. 



we can use the expression Eq. (|4"T|) to write the expecta- 
tion value of the same operator in the strip as 



(OH) 



I— I 

1 dz 1 



-Xo 



(O(z)) 



A° B \2- 



-X 



(69) 



Using \z\ = e -&( w+w "> 
ef u sin (jw) , we obtain 



(0(w)) 



and x = Im 



-X 



(70) 



where v = Im(u>) is physically the distance from the left 
boundary (v = 0) in the finite system (see the right-hand 
side of Fig. [9]) . Since the zero-temperature finite system is 
invariant under translations in it, the static ground-state 
expectation value of the local operator Jg(x)jf (x), in 
the finite system of length £ obtained from the transfor- 
mation Eq. ([67]) . is expected to behave as 



. 7T \ .7T 

2an( 7 *J/ 7 



(71) 



with x the distance from aboundary. 

A simple check of the above expression is that it re- 
produces, as expected, the correct power-law expectation 
value of the semi-infinite system Eq. (|66[) in the limit 
£ — > 00. By combining the above expression with the 
result of Eq. (|64|). we obtain 



G a0 = (J%(x)J p L (x)) GS 



sin 



(r) 



X < 72 > 



Notice that as long as the correlation function in the 
above expression is calculated for large enough x, i.e., 
larger than the microscopic length scales and the re- 
quired healing length for the contributions of the irrele- 
vant boundary operators to become negligible, the value 
of G a p should be independent of x. Practically, one can 
extract G a p by fitting the calculated ( Jr(x) jf (x))gs to 
the form of the universal sine function. Notice that the 



assumption that all wires have the same Luttinger pa- 
rameter g was not essential for the derivation of Eq. (|72p . 
This is because the generic form of Eq. ([fTl]) holds even 
if the wires a and /3 have different Luttinger parameters. 
The overall coefficient of ^-j^p however will depend 
on the two Luttinger parameters as well as the boundary 
condition B. 



The Eq. ([72]) is one of the key results of this paper. It is 
significant both from a fundamental and a practical point 
of view. From a fundamental viewpoint, we think about 
conductance when we have a state with currents flowing 
through the system from a source and into a drain. It is 
remarkable that one can construct a closed finite system, 
the ground state of which fully encodes the conductance, 
a quantity defined in open systems. From the practical 
point of view, the above expression is the basis for a 
method of calculating the universal conductance that we 
discuss in the next section where we show how to apply 
this continuum result to a lattice computation and how 
to implement the boundary conditions. 



VI. A METHOD FOR CALCULATING THE 
CONDUCTANCE 

In the previous section, we derived a general expres- 
sion that related the off-diagonal elements G a p of the 
conductance tensor to the static correlation functions 
(J^(x)J^(x)). The key to this relation is that both the 
conductance and the static correlation function above are 
proportional to only one coefficient that depends on the 
conformally invariant boundary condition at the junc- 
tion. This relationship can be used as the basis of a 
method for calculating the conductance. 

If we measure, numerically, the static correlation func- 
tion (Jr(x)Jl{x)) in a finite system obtained from the 
conformal transformation Eq. (foT)) . we can then fit the 
correlation function (a function of x) to the universal 
form Eq. (fTTj) and obtain the coefficient A°g ' . It is impor- 
tant to note that the universal expression is the asymp- 
totic form of the correlation function and the fit should 
be done in a region far away from the boundary. We 
need to consider large enough x so that the microscopic 
details become unimportant and the corrections due to 
irrelevant operators become negligible. 

To carry out this scheme, we need to compute 
(J^(x)J^(x)) in the ground state. The computation is 
amenable to the time-independent DMRG method. In 
the quasi-lD system obtained by having all wires run- 
ning parallel to one another, this correlation function is 
actually the expectation value of a local operator, which 
makes it a particularly easy quantity to calculate with 
DMRG. We need to answer two questions, however, be- 
fore we can proceed. 

1. Since we cannot model chiral creation and annihi- 
lation operators on a lattice, how can we model the 



operator J^{x) J^{x)\ 
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2. What exactly is the finite system that we need 
to put into a computer? In other words, given 
a tight-binding Hamiltonian for the original semi- 
infinite system, what is the Hamiltonian of the 
finite system for which the correlation function 
(Jff(x)J^(x)) has the universal form we derived in 
the continuum? 

To answer the first question, we make use of the follow- 
ing continuum relationship between the chiral currents 
and the total current and the total charge density: 



J" 



v{J%-Jt), N a = J% + J° L 



(73) 



Notice that since with electron-electron interactions these 
chiral currents are not the bare left and right-moving 
currents, the velocity v is not the Fermi velocity but the 
rcnormalizcd velocity of the charge carriers that can be 
calculated from the Bethe ansata^ and is given by the 
following expression at half-filling: 



yi-(v/2) a 

arccos (V/2) ' 



(74) 



where we have set the hopping amplitude to unity and V 
is the interaction strength in the bulk lattice Hamiltonian 
Eq. (|27p. From the Bethe ansatz, we can also obtain the 
effective Luttinger parameter in terms of the microscopic 
lattice Hamiltonian and, at half-filling, we have 

(75) 



2 arccos (-V/2)' 

Can anything go wrong if we apply the continuum rela- 
tion ([73")) to a lattice computation? The relation Eq. ([73")) 
is only exact in the continuum. However, we find that 
it gives very accurate results on the lattice for correla- 
tion functions between different wires a ^ (5. We will see 
this explicitly in the numerical studies of the following 
sections. The simple lattice calculations we performed in 
Sec. IIIII also shed some light on this issue. 

These off-diagonal correlation functions are those that 
we actually need for the conductance calculation, as in 
our approach, we first find the off-diagonal elements of 
the conductance tensor, which are simply proportional 
to A'jf . The diagonal ones can be later deduced from 
Eq. (|55|) . Now, we would like to find the off-diagonal 
elements of the conductance tensor G a p, and we use 
Eq. ([73")) as if it were an operator identity. Using Eq. ([75} , 
we can then write, for a ^ /3, 



(J a (x)J?(x)) = -v 2 ((J£(x)4(x)) + {J%{x)Ji{x))) , 
{N a (x)J^(x)) = v ({J£(x)4(x)) - {J%(x)jfc)) 



The correlation function (J^(x)J^(x)) is then simply 
given by 



1 



2v a V/3 



(J«( x )jP( x ))- — (N a (x)jP(x)), 



2vp 



where v a (vp) is the velocity in wire a (/?). Throughout 
this paper, we work with wires with the same Luttinger 
parameter and therefore v a — vp — v but the more gen- 
eral expression above can be used in cases where the wires 
have different Luttinger parameters. 

In writing the above expressions, we have made use of 
the fact that the chiral correlation functions of the form 
{Jr(%)Jb(x)) vanish in this finite system for We 
argued in section [V] that these correlation functions van- 
ish in the semi- infinite system (upper half-plane). Since 
the finite system corresponds to the strip obtained from 
the upper half-plane by a conformal transformation, they 
also vanish in the finite system. 

Note that if we have time-reversal symmetry, as in the 
case of the simple two- wire problem studied in section ITTll 
the second term in the above equation vanishes. This 
is because (nonchiral) currents J are odd under time- 
reversal and the densities N are even. If the structure 
at the junction does not include magnetic fluxes, time- 
reversal symmetry remains unbroken and we only need 
to calculate one current-current correlation function. 

The operators J and N are easily modeled in terms 
of the lattice creation and annihilation operators. Two 
important observations are in order. First, since the cur- 
rent operator is defined for a bond, we have to write an 
effective bond density so the current and density are as- 
sociated with the same position x. This is done by taking 
the average of the two site densities. The second obser- 
vation is that bosonization is done for the charge density 
fluctuations and the background charge should be sub- 
tracted to construct the operator N. With these two 
observations, we can explicitly write 

J (m + — ) = i{. c a ,m+l C ot,m ~ c tx,m c ct,m+l) > 

N a (m + -) = - (ra Q , m + n a . m+1 - (n a . m ) - (n a , m+ i)) . 

We now address the question of implementing the 
boundary conditions. To measure the expectation values 
above, we need the ground state of a finite system. What 
we need to do is write a lattice Hamiltonian that corre- 
sponds to the continuum system we obtained from a con- 
formal mapping. As discussed in the Introduction, the 
starting point is the Hamiltonian of the junction, which 
generically can be written as 



H = H] 



boundary 



bulk- 



The finite system has two boundaries on the left and on 
the right and as shown in Fig. 1101 generically, a Hamil- 
tonian 



H = Hn + Hi 



bulk 



(77) 



(76) 



Since the bulk theory is conformally invariant, we can 
write -Hbuik by using exactly the same terms as in i^buik> 
but just a finite number of them. The more subtle issue is 
to identify the boundary Hamiltonians Hq and Hi that 
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J / 



ndary 



H' = H Q + H! 



FIG. 10: (Color online) Given the bulk and boundary Hamil- 
tonians of a semi-infinite system, the goal is to obtain the two 
boundary Hamiltonians Ho and He of a finite system corre- 
sponding to the conformal transformation to the strip. 



MO) = So MO)- (78) 
As evident from Fig. [IT] the left and the right movers 

*out(0) = *h(0) . %n(£) = 



4 



*in(0) =*i(0) 



*outW = *LW 



give the correct conformally invariant boundary condi- 
tions. There are some details regarding persistent cur- 
rents and the parity of the number of electrons that we 
discuss later, but for now let us give a generic argument 
regarding the boundary conditions. 

The semi-infinite system that corresponds to the 
upper-half plane has one boundary and a conformally 
invariant boundary condition on this boundary, i.e., the 
real axis. Now, a conformally invariant boundary con- 
dition, by definition, does not change under a conformal 
transformation. The strip on the w = u + iv plane, how- 
ever, has two boundaries at v = and v = £. Let us 

review the conformal transformation w = — Inz. Writ- 
er 

ing z in polar coordinates as z = re % ^ gives v = -(f), 
so the left boundary comes from the positive real axis 
(0 = 0) and the right boundary from the negative real 
axis (cf) = 7r). Therefore, both these boundaries are ex- 
pected to have the same boundary condition as the real 
axis in the semi-infinite system. 

Now, in the semi-infinite system, this unknown confor- 
mally invariant boundary condition is stabilized by the 
structure and interactions of the junction. That is to 
say a cap with Hamiltonian i/boundary gives rise to the 
boundary condition B at x = 0. It is then a completely 
natural assumption that the same cap should give the 
boundary condition B at the left boundary of the finite 
system. This simply yields 

H-Q -^boundary- 

The boundary Hamiltonian on the other side (Hi), 
however, is generically different than Ho. Let us for con- 
venience represent the distance from the left boundary by 
x in the finite system as well as in the semi-infinite one. 
As argued above, we wish to have the same boundary 
condition on both sides. This does not mean, however, 
that we should place the same cap on both ends of the 
finite system. To understand the issue, let us consider 
the case of noninteracting electrons, where the boundary 
condition can be simply written in the form of a single- 
particle scattering matrix. 

Suppose the S matrix of the left boundary (x = 0) is 
So- As seen in Fig. [Til the incoming scattering states 
at x — are the right-moving fcrmions ipn(0) and the 
out-going states are the left- moving fcrmions ^l(O). The 
boundary condition can then be written in terms of So 



FIG. 11: (Color online) A generic finite system with two 
boundaries. For noninteracting fermions, the boundary con- 
ditions can be written in the form of scattering matrices that 
relate the in-coming and out-going states. The correspon- 
dence of these states with the chiral fermions is illustrated at 
both boundaries. 

switch roles at the two boundaries. Namely, the left 
movers (right movers) are the out-going (in-coming) scat- 
tering states for the right boundary and the in-coming 
(out-going) scattering states for the left boundary. If the 
scattering matrix corresponding to the Hamiltonian of 
the right cap Hg is Si, we have by definition 



Mi) = s t M*)- 



(79) 



However, we want to have the same boundary condi- 
tion at the two boundaries. This means that to write the 
boundary condition of the right boundary at x — I, we 
should write exactly the same expression as Eq. (|78[l but 
replace the argument x = with x = i, which yields 



MQ 



Comparing Eqs. ([79]) and 
important equation 



So 1 Mi) = 4 Mi)- (80) 

80| leads to the simple but 



s 



(81) 



Still working with a noninteracting system with the 
boundary condition (|78[) for the semi-infinite system, let 
us derive the boundary conditions of the finite system 
more carefully. The chiral fermionic creation and anni- 
hilation operators are in fact primary fields with scaling 

dimension — and transform as 
2 



' dw 

= ( Tz 



' dw 

*«H = ( - 



(82) 



By using the derivatives 

dw 
dz 
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we find that, at the two endpoints of the finite system, 

-e-f", (83) 

T 

— e-f". (84) 
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The minus sign above indicates that there could be 
subtle issues with the assumption of having the same 
boundary condition at both ends, which lead us to write 
Eq. Of course, notice that the minus-sign differ- 

ence between the transformed chiral fermions at the two 
cndpoints of the finite system is not inconsistent with 
Eq. ([5TH because we also need to take the square root 
of 4™ and H as seen from Eq. (|S2J. At this level, it 
is not clear, gencrically, which branch of the square root 
we should pick, but the two choices lead to the follow- 
ing two possibilities for the scattering matrix at the right 
endpoint: 

I : Si = Sq , II : Si = —S . 

The above relations can be thought of as two micro- 
scopic BCs, which in conjunction with other microscopic 
details such as the number of electrons in the system 
give rise to the correct BC in the continuum. Although 
we can microscopically implement both of these scatter- 
ing matrices, we work with the relation I above in this 
paper. Our approach is to use the relation I and choose 
the number of electrons such that, as discussed below, 
no persistent currents are generated in the finite system 
with this microscopic BC. 

The presence of persistent currents can give correc- 
tions to the CFT predictions. The reason persistent cur- 
rents need to be avoided is that the best region for fitting 
the correlation function (J^(x)J^(x)} is away from the 
boundaries and close to the center of the finite system 
where the effect of the irrelevant boundary operators is 
negligibly small. In the center of the finite system, the 

correlation function will be 0(—^) and since persistent 
1 

currents go as — , the corrections to (J^(x)J^(x)) , which 
are not accounted for in our formulation, become of the 
same — order as the value of the (J^(x)J^(x)) correla- 
tion function. 

As an example for how to avoid persistent currents con- 
sider a two- wire system with the boundary condition I, 
i.e., the system shown in Fig. |4]of Sec. IIIII We find that 
the persistent currents are avoided by having an even 
number of electrons. This can be seen by noting that for 
hopping amplitude at the junction t = 1 (a simple loop 
with a 7r flux and no actual impurity), the allowed mo- 
mentum levels exclude k = and are symmetric around 
it. Therefore, filling in an odd number of electrons leads 
to degeneracies. 

Having gained some intuition regarding the bound- 
ary conditions in the absence of electron-electron inter- 
actions, we now tackle the problem of constructing the 
Hamiltonian Hi of the right cap from that of the left, 
i.e., Hq — -ffboundary- We do this in two steps described 
below. 

First, consider the case where the interactions are only 
in the wires and -^boundary is quadratic. In this case we 
imagine turning off the bulk interactions. We then en- 
gineer a Hamiltonian Hi in a completely noninteracting 



system, which gives the scattering matrix Si = Sq at the 
desired filling. We now have a microscopic system with 
the same boundary conditions at both endpoints. Turn- 
ing on the interactions causes the system to flow away 
from the original noninteracting fixed point. The bound- 
ary conditions at both endpoints also flow away. How- 
ever, since they were the same initially, we expect them 
to flow together and reach a new BC. The two bound- 
aries are then expected to have the same BC even with 
interactions. 

The next step is to find symmetry transformations that 
automatically generate the desired Hamiltonian Hi from 
the given Ho. In what follows we argue that at half- 
filling, a combination of time-reversal and particle-hole 
transformations does the job. Having the boundary con- 
ditions implemented with such simple transformations 
strongly suggests that these transformation should also 
implement the correct BCs even when we have interac- 
tions in the junction itself as well as in the bulk of the 
wires. 

Our starting point for this procedure is an explicit ex- 
pression for the scattering matrix from Sec. IIII| which 
we repeat here for convenience: 

S(k) = -(T + e- lk )- 1 (T + e tk ), (85) 

where T is an M x M Hermitian matrix defining the 
boundary Hamiltonian as 

-^boundary = *o T ^o, 

with the notation defined in Sec. IIIII right above Eq. ([7)l. 
At half-filling, we have kp = — and the expression for 
the scattering matrix at the Fermi level reduces to 

Sf = -^4 • (86) 
I — i 

Note that it is important that the correct boundary con- 
dition is implemented at the Fermi level. Neither the 
conductance nor the behavior of the static (J^(x)J^(x)} 
correlation functions are strongly affected by what hap- 
pens deep inside the Fermi sea. It is evident from the 
form of the scattering matrix in Eq. ([5o]) that we need to 
change T — > — T to invert the scattering matrix. In other 
words, if we have 

Ho = *Sr * , Hi = tfjl^, (87) 

and we require Spi{kp) — S\{kp) at the half-filling, we 
must have 

r f = -r . 

This is one example where we in fact engineered the 
Hamiltonian Hi to implement the boundary condition I. 
It is very helpful to find symmetry transformations that 
automatically reproduce Hi from Hq. We argue that the 
transformation T — > — T for the quadratic Hamiltonian 
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Hq is a combination of time-reversal and particle-hole 
transformations. 

Let us call these transformations K and C, respec- 
tively. Since here we are considering spinless fermions, 
the time-reversal transformation is simply equal to the 
complex conjugation. Note that generically the time- 
reversal operator is anti-unitary and can be written as 
the product of the complex conjugation and a unitary op- 
erator. The particle-hole transformation simply switches 
the role of creation and annihilation operators. We can 
then write 

K(i) = -i, C{c) = c\ C*(c t ) = c. (88) 

Acting with the particle-hole transformation on the left 
boundary Hamiltonian of the form Eq. (1871) yields 

M 

c(H ) = c(*tr * )=5^c(r^4 )0Cfll0 ) 

= - ^ c a ,o c 0,O + const. 

Neglecting an unimportant constant, we can write the 
above transformation in a more compact form as 

c(h l ) = -*Sr?* = -^5*0, 

where we have used the Hermiticity of Tq in the second 
step. We can then write 

K(C(H a )) = -H = H e . (89) 

So far we have verified that for a simple junction de- 
scribed by the boundary Hamiltonian Eq. (|87j) . acting 
with K and C transformation on the boundary Hamil- 
tonian inverts the Fermi-level scattering matrix at half- 
filling. We conjecture here that this is a more generic 
statement, and acting with these transformations on 
the boundary Hamiltonian gives a microscopic boundary 
Hamiltonian that correctly implements the same bound- 
ary conditions at both endpoints. To support this conjec- 
ture, we test it on a more complex junction constructed 
by adding one more column of sites to the system. This 
calculation is presented in Appendix [A] 

The method developed in this paper can also be ex- 
tended to study systems with spin-i electrons. If there is 
no electron-electron interaction in the bulk of the wires, 
we can visualize the spin-up and down sectors as two 
distinct wires of spinless fermions. This construction 
can be used to engineer the Hamiltonian of the mirror- 
image junction from the scattering approach. In the pres- 
ence of interactions in the bulk, the Hamiltonian of each 
wire can be written as a sum of a charge-sector Lut- 
tinger liquid with velocity v c and Luttinger parameter 
g c and a spin-sector Luttinger liquid with velocity v s 
and Luttinger parameter g s . The conductance is then 
related to the the chiral current-current correlation func- 
tion (Jcr( x )Jcl( x )) m * ne charge sector. To measure this 



correlation function, we need to calculate v c in terms of 
the microscopic parameters such as the strength of the 
on-site Hubbard interaction Un^n^, where f (4-) indicates 
spin up (down)i^ This is a bulk property and can be 
calculated numerically or via the Bethe ansatz. The cor- 
relation function {Jcr{x)J^ l {x)) can then be computed 
numerically using the following relations: 

Jj? + J? = v c (J« R - N? + N? = J« R + J? L , 

which are analogous to Eq. (|73p . These relations lead to 
a similar equation to Eq. (|76|) . 

One final comment is in order regarding the efficient 
computation of the required correlation functions with 
DMRG. When interactions are confined to the junction 
itself, i.e., in the absence of interactions in the bulk of the 
wires, it is in some cases possible to perform a canonical 
transformation, which decouples some of the degrees of 
freedom in the system as recently shown in Ref. l59l . Such 
transformations when possible can reduce the dimension 
of the local Hilbert space in DMRG calculations and lead 
to significant performance improvements. 

VII. NUMERICAL TESTS AND BENCHMARKS 

A. Noninteracting Y junction 

The first set of numerical tests and benchmarks we per- 
form is with a noninteracting Y junction. The system we 
consider has three noninteracting fcrmionic chains with 
hopping amplitude of unity. At the junction, each chain 
is coupled to the other two chains with hopping ampli- 
tude t and the loop formed at the junction is threaded 
with a magnetic flux $. The gauge-invariant flux $ can 
be modeled in many different ways in the tight-binding 
Hamiltonian such as, for example, the following bound- 
ary Hamiltonian: 

/ te^ A 
H = te- J * t * - (90) 

\ t to) 

Using Eq. © and the Landauer formula ©, we obtain 
the following exact expression for the conductance of this 
simple junction: 

G 1221 = 4f(l + f±2tsin$) 
12,21 1 + 6t 2 + 9t 4 + At e cos 2 $h~ { ' 

We can now use exact diagonalization of the single- 
particle Hamiltonian, which can be written as an 37V x 3 A 
matrix (with N the number of sites in each leg and 3A 
the total number of sites in the system), to obtain the 
(J a (x)jP(x)) and (N a (x)J^(x)} ground-state correla- 
tion functions. Notice that since the total Hamiltonian is 
quadratic, Wick's theorem [Eq. (|15p] holds. Let us label 
all the 3A sites in the multi-wire system in some order 
as, for example, in Fig. I12[ i.e., assign an annihilation 
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operator c a , a — 1 . . . 37V to every site. All we need to 
calculate are then fermionic Green's functions of the form 

C(a,(3) = (cicp). (92) 



a = 4 




- 1 



FIG. 12: A labeling of system sites for exact diagonalization. 

Diagonalizing the 3N x 3 AT single-particle Hamiltonian 
gives single-particle wave functions |* e ) = J2 a 4> e a \a) 
where \a) is a basis state of the single-particle Hilbcrt 
space with site a occupied and all other sites empty and 
e represents the energy of the single-particle level. These 
wave functions are obtained directly from exact diagonal- 
ization. At a given filling is, the many-body wave function 
is a Slater determinant of 3vN low-energy single-particle 
wave functions. The wave functions of the filled levels 
are then <jf a for e < ep where ep is the Fermi energy. We 
can then write 

C(«,jJ)=E^. (93) 

e<e F 

We calculate, through exact diagonalization, the cor- 
relation function (J^(x) J^(x)) as explained above and 
plot their logarithms versus In [— sin for different 

values of the hopping t and for two different fluxes $ = 
and <& = 7r/2. The results for $ = are shown in Fig. [13] 

As expected from the analysis of previous sections, far 
away from the junction, these plots are lines of slope —2. 
The conductance can be extracted from fitting the data 
in Fig. [T3] and, as shown in Fig. [TJ] it agrees with the 
exact conductance Eq. (|91[) with very good accuracy. 

B. Junction of two interacting quantum wires 

The next benchmark we perform is for two interacting 
wires with an impurity. As shown in Ref. this system 
exhibits universal behavior. For attractive interactions, 

any impurity should heal at the RG fixed point, leading 

2 

to a conductance of Q%r, while for repulsive interactions, 
any impurity should effectively cut the chain resulting in 
zero conductance. 

First, we show benchmarks for the attractive case. 
Here, of course, we are dealing with a strongly correlated 
system and exact diagonalization of the single-particle 
Hamiltonian can not be used. Due to the exponen- 
tially large dimension of the Hilbert space, we cannot 
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FIG. 13: (Color online) The inset shows the noninteracting Y 
junction used for the numerical test of the method. The thick 
bonds at the junction have a different hopping amplitude t 
than the thin bonds in the bulk of the wires, the hopping 
amplitude of which is 1. The static correlation function is 
for a finite system of length I — 99 (100 sites in each wire) 
for $ = and five different values of the hopping amplitude 
t. The lines drawn through the data points show the CFT 
prediction with the exact conductance Eq. (|91|) . The strong 
agreement with the numerical data confirms that our method, 
i.e., extracting the conductance from fitting the correlators, 
yields accurate results. 

1 , , , , , , , — I 



0.8 




— Gl2, 


$ = 


0, Exact 






VG12, 


* = 


0, Numerical 


0.6 




— Gl2, 


4> = 


5, Exact 


Oo.4 




□Gl 2 , 


<I> = 


Numerical 






-~G ai , 


* = 


|, Exact 


0.2 




oG 21 , 


<1> = 


7?, Numerical 



qI 1 1 1 a^r 1 1 1 * — i 

-3 -2 -1 1 2 3 4 5 

t 

FIG. 14: (Color online) The conductance of the noninteract- 
ing Y junction calculated with our method versus the exact 
solution. 

exactly diagonalize the many-body Hamiltonian cither. 
The DMRG method is an ideal tool for performing these 
calculations, and the results presented in this section are 
obtained by this method. 

The results of these calculations are shown in Figs. [T51 
and 1161 for the Luttinger parameter g = 2.0 and g = 6.0, 
respectively. For the larger Luttinger parameter g = 6.0, 
we are very close to the phase separation point V = 2 and 
we need a larger healing length. The universal behavior 
is evident from the fact that the value of the hopping 
amplitude at the junction has no effect on the asymptotic 
form of the correlation functions. 

Let us now consider the case of repulsive interactions. 
This is a special case because the theoretical prediction 
for the conductance is G12 = so we do not expect to 
be able to fit the correlation function to the CFT form. 
The behavior of the correlation function, however, sheds 
light on the nature of the problem. As noted earlier, the 
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FIG. 15: (Color online) The correlation functions are plotted 
for two values of hopping at the impurity t = 0.8 and t — 0.6 
for the Luttinger parameter g = 2.0. The strong agreement 

2 

with the CFT prediction of G — g^ supports the correctness 
of our method. 
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FIG. 17: (Color online) The correlation functions are plotted 
for hopping at the impurity t — 0.3 and the Luttinger param- 
eter g — 0.65. Since the universal conductance vanishes, we 
only observe the nonuniversal subleading corrections. 
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an open problem for three interacting wires. 




FIG. 16: (Color online) The correlation functions are plotted 
for two values of hopping at the impurity t = 0.8 and t — 0.6 
for the Luttinger parameter g = 6.0. Here a larger healing 
length of I = 90 is needed to observe the universal behavior 
and verify the CFT prediction. 



CFT prediction for the universal form of the (J^Jf) is an 
asymptotic form, and we expect subleading corrections 
coming from irrelevant boundary operators. 

In the repulsive case, since the conductance is zero, 
(JrJl) wm decay faster than [sin(xj)] . As seen in 
Fig. [T7] for g = 0.65 < 1, where the coefficient in front of 
the leading term is zero, we only observe these subleading 
corrections. As a function of In [sin(xj)], In (Jjjj£) de- 
cays with a larger negative slope than —2 and in addition 
exhibits a 2kp oscillatory behavior. 

The studies of this section verify that the method de- 
veloped in this paper correctly reproduces established re- 
sults in the presence or absence of electron-electron in- 
teractions. In the next section we apply the method to 



VIII. UNIVERSAL CONDUCTANCE IN AN 
INTERACTING Y JUNCTION OF QUANTUM 
WIRES 

In this section, we employ the method developed in this 
paper to study a Y junction of interacting quantum wires 
with spinless electrons, each connected to a distinct site 
on a loop at the junction as seen in Fig. [3] The tunneling 
amplitude between the wires is t and the flux through 
the l oop is d>. This system was analytically studied in 
Refs.ka andl46l where several nontrivial fixed points were 
identified. A similar anal ysis for a Y junction with spinful 
electrons is given in Ref. roll 

In this work, we focus on the Luttinger parameter 
range 1 < g < 3. The RG flow diagram for this range, 
which was conjectured in Ref. l46l , is shown in the right- 
hand side of Fig. [18] We have two types of fixed points: 
the chiral fixed points \+ an d X- in the presence of a 
time-reversal-symmetry-breaking flux and the M fixed 
point for unbroken time-reversal symmetry. 



A. Chiral fixed point 

For the chiral fixed point, the time- reversal symmetry 
is broken by a magnetic flux and therefore G\i ^ G2i- 
It was predicted in Ref. [46| that these conductances ac- 
tually have a different sign and depend on the Luttinger 
parameter as follows: 
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FIG. 18: The conjectured RG flow diagram for the Y junction 
of Fig. [31 given in Ref. [4(| for the Luttinger parameter in the 
range 1 < g < 3. 



FIG. 19: (Color online) A plot of the ln{ J^x) Jl(x)) for 
g = 1.5 and g — 2.0 and two different hopping amplitudes. We 
observe the universality of the conductance and good agree- 
ment with the theoretical predictions. 



This conductance is universal, i.e., it is independent of 
the hopping amplitude t at the junction for nonvanish- 
ing t. It is worthwhile to emphasize that the emer- 
gence of universal conductances strictly relies on the 
presence of interactions, i.e., by taking the noninteract- 
ing limit g — > 1, we shall not expect to recover correct 
results for noninteracting cases £2- In the RG language, 
the noninteracting point is a marginal fixed point and 
each value of hopping amplitude t gives a unique conduc- 
tance as in Eq. (|91|) . Because we are studying junctions 
with Z3 symmetry, the conductances obey the relations, 
G\2 = G23 = G31 and G 2 i = G 32 = G13. In the time- 
reversal symmetric case (the M fixed point), the conduc- 
tances in addition follow = Gji. To verify these pre- 
dictions, we numerically calculated the correlation func- 
tion, (jj i (x)Jl(x)}, for two values of t and two values of 
g in a system size of I = 40. As shown in Fig. (THl the 
numerical results strongly agree with the prediction for 
G\2- The consistency serves both as the first numerical 
verification of the prediction Eq. (|M)) as well as a highly 
nontrivial benchmark for the correctness of our method. 

It is quite remarkable that we can verify the G12 con- 
ductance in a system as small as I = 40. Note that due to 
the asymptotic form of the relation between the univer- 
sal conductance and the computed correlation functions, 
the results are reliable once there is convergence in the 
system size. Also, one needs to check for the conver- 
gence of the DMRG method by increasing the number of 
states used in the calculations to make sure the DMRG 
results are in fact quasi-exact. We show some of these 
convergence tests in Fig. [201 

In the last part of this section, we present the numeri- 
cal verification of the G21 conductance. It turns out that 
for G21, a longer healing length is needed and we are 
able to see the agreement with the theoretical prediction 
in a system of I = 60. It is remarkable that the correla- 
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FIG. 20: (Color online) In the left panel, we check the conver- 
gence of the correlation function with increasing system size 
for g = 1.5 and t = 0.7. As seen, the correlation functions 
collapse for £ — 40 and I = 60. In the right panel, we check 
the DMRG convergence as the truncation error is reduced by 
increasing m. The data are for g — 2.0 and t — 0.7 in a 
system size of £ — 60. 

tion function here actually changes sign and the universal 
long-distance behavior, corresponding to a positive con- 
ductance, has opposite sign to the nonuniversal values of 
the correlation function near the junction. 

Due to the sign change of the correlation function, here 
we plot the logarithm of minus the correlation function in 
a region far away from the junction (where the argument 
of the logarithm is positive). The results are shown in 
Fig.[2T1and. once again, there is good agreement between 
the data and the theoretical prediction. 



B. M fixed point 

Let us now turn to the conductance of the time-reversal 
symmetric M fixed point for which theoretical predictions 
do not exist. The benchmarks in the previous section and 
the study of the chiral fixed point for three interacting 
wires increases our confidence in the correctness of the 
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FIG. 21: (Color online) The universal G21 conductance for 
g = 2.0 and two different values of the hopping amplitude 
t — 1.0 and t — 0.7. Notice the sign of the argument of the 
logarithmic function labeling the vertical axis. Here, G12 is 
positive and the correlation function (J^ l (x)jj j (x)} actually 
changes sign at a certain distance from the boundary. 



method developed here. We are now ready to apply the 
method to a thus-far unsolved quantum impurity prob- 
lem. 

We perform the calculations of the M fixed-point con- 
ductance for three values of the Luttinger parameter g 
and two values of the hopping amplitude t. The indepen- 
dence from the hopping amplitude can be seen at t = 60. 
The results are shown in Figs.[22l[23l andl24l respectively 
for g = 1.5,2.0,2.5. 
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FIG. 22: (Color online) The universal G12 = G21 conductance 
for g — 1.5 and two different values of the hopping amplitude 
t = 1.0 and t = 0.7. 

To obtain the conductance, we fit the data for a slightly 
larger system of I = 70. The three correlation functions 
are plotted in Fig. [51] in addition to the best fit in the 
asymptotic region for each value of g. We obtain the 
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FIG. 23: (Color online) The universal G12 = G21 conductance 
for g = 2.0 and two different values of the hopping amplitude 
t = 1.0 and t = 0.7. 
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FIG. 24: (Color online) The universal G12 = G21 conductance 
for g = 2.5 and two different values of the hopping amplitude 
t = 1.0 and t = 0.7. 

following conductances for the M fixed point: 

G 12 (g = 1.5) = G 2 i(g = 1.5) « -0.55^-, 

h 

G l2 (g = 2.0) = G 21 (g = 2.0) « -0.62^-, 

h 

G 12 (.g = 2.5) = G 21 (g = 2.5) w -0.665 ^-. 

Is there a pattern to these values of the conductance? 
It is well known that when a Luttinger liquid with at- 
tractive interactions (g > 1) is coupled to Fermi- liquid 
leads, the ge 2 /h conductance renormalizcs to e 2 //t i 63 i 64 
One way to think about this phenomenon is to assume 
that we have an effective contact resistance at the inter- 
face between each interacting wire and the noninteracting 
leads. If this contact resistance is l/g c , we must have 

l/g + 2/g c = l g c = 2/(l-g- 1 ). 
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FIG. 25: (Color online) Fitting the correlation function to the 
CFT prediction to extract the conductance. The data in this 
figure are for t = 0.7 and £ = 70, with three different values 
of g. 



For the multi-wire conductance tensor, we can consider 
the following transformation, which, if the effective con- 
tact resistance is the correct mechanism, renormalizcs 
the conductance in the presence of nonintcracting leads. 
Here, we think of the transformation 

G = (1 + G/^r 1 G 

as a mathematical mapping of the conductance tensor, 
which can be inverted as 

G = G(l-G/. 9c ) _1 . 

It turns out that for all the other fixed points of the Y 
junction as well as in the two- wire case, G is independent 
of the Luttinger parameter^ We conjecture here that 
this is also true for the M fixed point. We can check 
explicitly that for the three values of g studied above, 



G 



with 7 w 0.42. This conjecture leads to a general form 
for the conductance of the M fixed point as a function of 
the Luttinger parameter and one dimensionlcss number 
7 w 0.42 as 




G12G?) = G 21 (g) 



257 



2.9 + 37 - 357 h 



As stated above, this conjecture works to a good ap- 
proximation for the three values of g studied above nu- 
merically. The functional form of the dependence of the 
conductance of the M fixed point on the Luttinger param- 
eter is plotted in Fig. [Ml This simple plot, which had 
remained elusive since the prediction of the existence of 
the M fixed point in^, is the main result of this section 



and an important achievement for the method developed 
in this paper. 

Notice that the measured value of 7 differs from 4/9 « 
0.444 by about 4.8%. There are several sources of er- 
ror in our evaluation of 7. The DMRG truncation error 
is very small in this case and the calculated correlation 
functions are quasi-exact. The error from fitting the data 
is around 2% and there are errors from finite-size correc- 
tions that with our system sizes are of the order a few 
percent. It seems plausible that the renormalized conduc- 
tance of the M fixed point, when attached to noninteract- 



ing leads, is equal to G12 



4„2 



/h as indicated by the 



results of Ref . |38| which were obtained by an approximate 
functional renormalization group scheme. This strongly 
suggests a simple understanding of the M fixed point con- 
ductance because §e 2 //i is the largest conductance IG12I 
we can have, with time reversal symmetry, for three non- 
interacting quantum wires connected by a unitary scat- 
tering matrix. The ^-dependence of the M fixed point 
conductance is such that when attached to noninteract- 
ing leads the system acts as an effective scattering matrix 
with the largest possible conductance G12. 



IX. CONCLUSIONS AND OUTLOOK 

In summary, we have developed a method in this pa- 
per that allows the computation of the conductance ten- 
sors of rather arbitrary quantum junctions with standard 
numerical techniques such as time- independent DMRG. 
We presented specific studies of an interacting Y junction 
and resolved an outstanding open question regarding the 
conductance of the M fixed point. 

Our generic method has only one practical limitation. 
Namely, the size I of this system must be large enough 
for the universal scaling behavior to emerge. Two factors 
affect this minimum length. First, an intrinsic healing 
length required to observe the continuum long-distance 
Luttinger-liquid behavior. For large attractive interac- 
tions (large g), this intrinsic healing length is large. Sec- 
ond, the scaling dimension of the irrelevant boundary 
operators. If these operators arc only slightly irrelevant, 
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we need to move away a larger distance from the bound- 
ary to observe the fixed-point behavior. The D p fixed 
point of the Y junction, for example, which occurs for 
3 < g < 9, presents these difficulties. The scaling dimen- 
sion of the leading irrelevant boundary operator increases 
with g, but making g larger brings the Luttinger liquid 
close to the phase separation point and increases the in- 
trinsic healing length. We have attempted to verify the 
theoretical prediction for the D p fixed point but have not 
been able to reach the required length with our compu- 
tational resources. As demonstrated in this paper, how- 
ever, this practical limitation is by no means a generic 
problem. 

The method developed in this paper opens up the pos- 
sibility of applying the powerful technology of DMRG 
and matrix product states, which in the last decades have 
been successfully applied to quasi-f D systems such as lad- 
ders, to transport calculations in the presence of electron- 
electron interactions. For decades, quantum transport 
was formulated for noninteracting electrons, and interac- 
tion effects were treated by mean-field and approximate 
methods. The results of this paper make it possible to 
utilize the quasi-exact numerical methods such as time- 
independent DMRG to study the transport properties of 
strongly correlated systems. 

Many interesting questions can now be answered with 
standard computations. A long-standing open question is 
the conductance of a chain of interacting spinful fermions 
at the nontrivial Kane and Fisher fixed point . 53 ' 65 ' 66 The 
method also makes it possible to study fixed points of 
more than three quantum wires. For two wires of spin- 
less fermions, there are only two interacting fixed points, 
while the Y junction presents a much richer flow dia- 
gram with multiple fixed points. What are the possible 
universal conductances with four, five wires, etc.? Are 
there hierarchies or other mathematical structures? 

Another direction of interest is to identify the uni- 
versality class that particular junctions, which are ex- 
perimentally accessible, fall into. Transport measure- 
ments are done on a multitude of organic molecules 
such as fullerenes, for example. Effective tight-binding 
models can be derived for such molecular structures via 
first-principles calculations. Our method combined with 
DMRG computations is then helpful for extracting the 
universal behavior of such molecular junctions. 

In summary, the technology developed in this paper 
provides a systematic way to study the effects of strong 
electron-electron interactions in the transport properties 
of quantum impurity problems and molecular electronic 
devices. In this framework, and with the help of efficient 
DMRG computations, strong correlations can be studied 
in a quasi-exact manner without resorting to mean-field 
or perturbativc treatments. 
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Appendix A: Symmetry transformation for a 
two-column noninteracting junction 

In this appendix, we test the Hr = K{C{Hl)) con- 
jecture of Eq. (j89]) for a two-column noninteracting junc- 
tion. The boundary Hamiltonian for the junction, which 
is shown in Fig. [271 is as follows: 

h l = *Jr * + *L!r_i*_i + v^T^o + *jTt*_ l5 

, .(Al) 

where the two T matrices arc Hcrmitian and T is an 
arbitrary M x M matrix. Let us first calculate the S 
matrix. We have 



(A2) 
(A3) 

1, we 



eip-! = r.iV'-i + Tip , 

eipo = r ip + r f v.'-i - ipi. 

We work at half-filling where e = 0. From Eq. 
then obtain 



tp-i = - (r_i) re- 
plugging this into Eq. (|A3|) and replacing the wave func- 
tions ip with scattering plane waves gives 

T (A out +A- m ) - (T^y'TiA^t + A^) 
- i(A out - A in ) = 

from which we get the following scattering matrix: 

S = 



i - r + t t (r_ 



T 



i + 



ro-Ttp.i) X T 



We have seen already that the transformation KC 
takes Tq to — Tq. It similarly changes r_i to — 
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In order for this transformation to work, we need to 
have Tt(r_i) _1 T -> -Tt(r_!) _1 T under KC. We 
will show this by arguing that the matrix T goes to — T. 
Let us begin by applying the particle-hole transformation 
to the T term tfljT^o + #jTt#_ x . We have 

C^liT^o + *o Tt *-i) 

= -J2{ T u ,c -u c h' +T h' c v c -ij') 

id' 

and, therefore, 

= ^LiC-Tj^o + ^C-rt)*-!. 

So, we find that indeed Tt (r_i) _1 T -> -T+ (r_i) _1 T. 



where is the boundary state on the real axis. 

The key to this derivation is calculating the expec- 
tation value of O on a semi-infinite cylinder of radius 
R seen in Fig. [25] in two different ways. First, we use 
a conformal mapping from the upper half-plane to the 
semi- infinite cylinder. Consider w = feff • It is easy to 
show that for z on the real axis, ww* = 1 and the upper 
half-plane is mapped to the outside of a unit disk. It is 
then easy to see that the conformal mapping 

R fl-iz 

w(z) = — In 

w 2tt \l + iz 

takes the upper-half plane to the semi-infinite cylinder. 
By applying this transformation to the expectation value 
of O on the semi-infinite plane, we obtain 



(0(y)) 



A° 



2tt\ 
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Xo 



eM-^X y/R) 
[1 - cxp(-47n//i?)] 
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Appendix B: Relationship between A^f and the 
boundary state \B) 



In this appendix, we review the notion of the highest- 
weight states and derive Eq. (|62[) . which explicitly re- 
lates the coefficient A°g to the boundary state and the 
highest-weight states of some primary operators in the 
CFT. The derivation in this appendix is merely a review 
of the generic result obtained in Ref. and is presented 
for completeness. 

In field theory, it is customary to quantize the theory 
for constant times. The Hamiltonian of the system is the 
generator of time translation. In CFT, however, because 
of scale invariance, it is convenient to quantize the the- 
ory on a fixed circle in the complex z = r + ix plane. 
This is called radial quantization and the generator of 
scale transformations is called the dilatation operator D, 
which can be written as the integral of the radial compo- 
nent of the stress tensor. 

The operator-state correspondence in CFT is formu- 
lated in the radial quantization framework. There is a 
vacuum state |0) in radial quantization, and acting by 
the operator corresponding to a scaling field on the vac- 
uum gives a state 

|0) = 0(O,O)|O). 

It can be shown that these states are eigenstates of the 
dilatation operator DJ^ The state constructed in this 
manner for a primary field is called the highest-weight 
state of that field. 

Let us now turn to the derivation of Eq. (|rJ2"j) follow- 
ing Reff££. We would like to show that if (0(x)) = 
Ag (2x)~ x ° for a scaling operator O living on the upper- 
half plane, we have 
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FIG. 28: The expectation value of O is calculated on the 
semi-infinite cylinder of radius R in two different ways. 

Now a second way to calculate this quantity is directly 
on the semi-cylinder and by using the boundary state. 
If the transfer matrix in the y direction (parallel to the 
cylinder axis) is exp(— H) for a Hamiltonian we can 
write 



(0(y)) 



lim 

T-HX 



(0\ e -( T -y) n O(y)e- yH \B) 
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\B) 



where |0) is the ground state of T~L. 

We now insert the resolution of the identity between 
0(y) and e~ v ' H in the above expression. Since O is a 
primary operator only states that are not annihilated by 
O enter the sum. The highest- weight state O\0) = \G) 
has the smallest eigenvalue — and, in the limit of large 
y, will be the leading term. We then obtain 



(0(y)) = 



(0\O\O)(O\B) 
(0\B) 



-£OV 
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where e = 2ttX /R~ We now argue that (0\O\O) be- 
haves as (2tt/R) x °. This can be seen from calculating 
(0|O(yi)O(?/2)|0) on a cylinder. By mapping the full com- 
plex plane to a cylinder of radius R. we obtain^! 



(0|O(yi)O(ifc)|0) = 



— sinh^(yi - y 2 ] 

7T R 



-2Xq 



In the limit of (y± —1/2) — > 00, the above correlation 
function reduces to 



(0\O( yi )O(y 2 )\0) 



R 
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Alternatively, we can write an expression for the above 
correlation function by inserting the resolution of the 
identity between Q(yi) and 0(3/2) and taking the limit 
of (yi — 2/2) — ^ °Oj which yields 

(010(2/00(2/2)10) ~ |(0|O|O>| 2 e- Eo ^-^). 



Comparing Eq. (|B1[) (in the limit of large R) and Eq. (|B2[) 
and making use of the above expression Eq. (|B3[) gives 
the following final relationship between the coefficient Ag 
in terms of the boundary state: 



(gig) 
(016) ' 



We then find by comparison that 

(0|O|O) ~ (27r/i?) Xo . (B3) 
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